General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



NASA Contractor Report 158944 

N 7 8- 334 79 

Unclas 
33826 

CREATION OF LUMPED PARAMETER THERMAL MODEL BY 
THE USE OF FINITE ELEMENTS 


(NASA-CR-158944) CREATION OF LUMPED 
PARAMETER THERMAL MODEL BY THE USE OF FINITE 
ELEMENTS (Sperry Support Services) 210 p HC 
A10/MF A01 CSCL 20K 

G3/39 


SPERRY SUPPORT SERVICES 
Huntsville, Alabama 35801 


NASA Contract NASI -14999 
October 1978 


NASA 

National Aeronautics and 
Space Administration 

Langley Research Center 

Hampton, Virginia 23665 



VOLUME I 


SECTION PAGE 


1.0 INTRODUCTION . . 1-1 

1.1 METHOD 1-1 

1.2 USER PROCEDURES . . . 1-3 

1.3 THERMAL NODE AND CONDUCTOR REQUIREMENTS 1-5 

1.4 MODELING TECHNIQUES 1-8 

2.0 THERMAL MESH GENERATION - SPAR PROGRAM 2-1 

2.1 CAPABILITY 2-1 

2.2 LIMITATIONS 2-3 

2.3 IN PUT/ OUTPUT DESCRIPTION . 2-3 

3.0 FINITE ELEMENT TO LUMPED PARAMETER CONVERSION - CINGEN 

PROGRAM . 3-1 

3.1 CAPABILITY AND METHOD 3-1 

3.2 LIMITATIONS 3-6 

3.3 INPUT/OUTPUT Description . . . 3-7 

4.0 INTERACTIVE GRAPHICS - GEOMPLT PROGRAM 4-1 

4.1 CAPABILITY 4-1 

4.2 LIMITATIONS . 4-14 

4.3 INPUT/OUTPUT DESCRIPTION 4-14 

5.0 SAMPLE PROBLEM 5-1 

5.1 SPAR MODELING 5-1 

5.2 CINGEN MODELING 5-1 

5.3 MODEL CHECK-OUT 5-1 

6.0 REFERENCES . . . 6-1 

APPENDIX A - SPAR USERS GUIDE A-l 

APPENDIX B - GEOMPLT USERS GUIDE B-1 


VOLUME 2 


APPENDIX C - CODE 


Section 1 


INTRODUCTION 

Thermal analysis by the lumped parameter finite difference technique is 
the most widely used and accepted approach to the solution of complex, 
mixed-mode, heat transfer problems. Computer programs such as CINDA^, SINDA^, 
and MITAS 3 are widely distributed and have found diverse application 
in aerospace and commercial industries. 

In the finite difference technique, the thermal network is represented 
by an analogous electrical network. The development of this network model, 
which is used to describe a physical system, often requires tedious and 
menial data preparation and checkout by the analyst. The data preparation 
and checkout can be greatly reduced through the use of the computer programs 
described in this report. These programs automatically develop the 
mathematical model and associated input data and graphically display the 
analytical model to facilitate model verification. 

Three separate programs are involved which are linked through common 
mass storage files and data card formats. These programs are SPAR 
(processors TAB, TELD and AUS), CINGEN and GEOMPLT. The objectives of the 
programs were to: 

1. Develop thermal models for the MITAS II thermal analyzer program 

2. Produce geometry plots of the thermal network 

3. Produce temperature distribution and time history plots. 

1.1 METHOD 

The basic approach used to generate a thermal model (MITAS data) was to 
develop an equivalent finite element model, convert this model 
to a lumped parameter thermal network, and display the network graphically 
for verification. Figure 1 illustrates the basic program flow diagram. 

A lumped parameter thermal model inherently lacks the geometric data 
necessary for graphical display. A finite element model, however, 
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Figure 1, Integrated Flow Diagram SPAR/ Cl NGEN/GEOMPLT 
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contains all of the geometric and physical data necessary to create a graphical 
displa,' . Furthermore, there are numerous finite element mesh generation 
programs which can be utilized to develop large analytical models with a 
minimum of data input and user preparation. Some of these programs are 
C00NS4, C00N3D5, ICES T0P0L0GY6, and SPAR TAB/ELD?. The latter program was 
selected for thermal model mesh generation since a modified version 
developed by Sperry has the capability of producing both NASTRAN^ and SPAR 
finite element thermal or structural models, and the simplicity of use was 
judged superior to the other programs. 

Having produced a finite element model using the SPAR formatted data 
cards, the second program (CINGEN) is used to convert the finite element 
model into a lumped parameter model. The output of this program is MITAS 
nodal and conductor data cards. 

Finally the resulting MITAS data and the SPAR finite element data are 

q 

processed by the GEOMPLT interactive graphics program for visually dis- 
playing the thermal conductor network and the finite element mesh. The 
analyst would normally view the model from several angles and "partial 
views." If model refinement or error corrections were required, the entire 
process may be repeated or the local file containing the data may be edited 
using the computer system text editor processor. When the nodal and 
conductor data are correct, they may be combined with the necessary MITAS 
control cards for processing in a thermal analysis. 

An option for post processing thermal data in the form of XY plots is also 
available in the GEOMPLT program. 

1.2 USER PROCEDURES 

Detailed input and output requirements for each program are described 
in Section 2, 3 and 4 of this report. The following discussion is an over- 
view of the procedure recommended for developing MITAS or SINDA data. 

a. When using the SPAR program for mesh generation, three processors 
must be executed. The TAB (CDC-SPARPUNCH) or TABX (UNIVAC) processor 
is used to generate joint (grid point) locations for the finite 
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element model. The user is required to set the "ONLINE = 3 or 4" flag 
which indicates that GRID cards will be punched. The TELD 
{CDC-SPARPUNCH or UNIVAC) processor is used to generate the 
NASTRAN finite element connection cards. The SPAR TELD (CDC or 
UNIVAC) processor can be used to obtain conduction, convection, 
and radiation elements. The user is required to RESET CINGEN=1 
to punch the element cards. 

b. Any other finite element mesh generation program may be used provided 
that standard NASTRAN formatted BULKDATA cards are punched or that 
the SPAR thermal element punch format is used. 

c. If convection or radiation is present in the thermal problem, the 
coefficients h c and eaF must be manually inserted into the finite 
element PUNCH file or CINGEN input deck. 

d. The user is now ready to execute the CINGEN program. The input 
cards are placed on TAPE4 (CDC) or file 4 (UNIVAC). This may be 
accomplished by renaming the PUNCH file resulting from the SPAR 
execution or loading the cards into a local (temporary) file 
(TAPE4). The CINGEN program will read this file and compute 
the lumped thermal characteristics of each finite element. This 
data will then be punched on cards in the appropriate MITAS (CDC) 
or SINDA (UNIVAC) format. In addition, a mass storage file 
TAPE29 (CDC) or 29 (UNIVAC) is written for use by the graphics 
program. 

e. When file TAPE4 and TAPE29 have been attached to a demand RUN using 
a Textronix terminal, the GEOMPLT program may be used to display 
either the finite element model or the thermal model or both. 

Labeling is provided for conductor and node type and number, partial 
views may be plotted and each view may be rotated to any viewing angle. 
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f. Finally, the user may create an output file from MITAS which contains 
tabular Node-Temperature or Time-Temperature data from the thermal 
analysis. This data can then be displayed by the GEQMPLT program 
as an X-Y plot. Node-Temperature plots will be scaled along the 
abscissa proportional to the geometric relationship of the nodes. 


1.3 THERMAL NODE AND CONDUCTOR REQUIREMENTS 


This section describes the type of lumped parameter thermal nodes 
and conductors which can be defined by SPAR finite elements and 
converted to MITAS data by the CINGEN program. The description 
includes the characteristic input for ^ach type. 


Diffusion Node - 


• Diffusion nodes are located at the geometric centroid of an 
K or F-type finite element. 

• K-type elements include bar (K21), triangular (K31) and 
quadrilaterial (K41) plates, wedge (K61) and hexahedron (K81) 
solids. F-type elements are restricted to bars (F21). 


• Positive values of density (p) required. 


Negative values of specific heat (C„) or conductivity (K) 
produce a reference to an ARRAY table. 


• ARRAY tables are not automatically generated and must be 
supplied by the user. 


Arithmetic Node - 

• Arithmetic nodes are located at the geometric centroid of any 
K or F-type finite element. 

. Any K-type element (K21 , K31 , K41 , K61 , or K81) or F-type (F21) 
may be used to generate an arithmetic node. 

• If area of K21 or thickness of K31 or K41 is input as 0, 
a surface arithmetic node is generated. 

• Negative value of density required. 

. Arithmetic surface nodes are generated for C and R - type finite 

elements provided these elements are not applied to a surface element. 
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Boundary Node - 


• Boundary nodes are located at the geometric centroid of any K or 
F-type finite element. 

. Any K-type element (K21 , K31 , K41 , K61, or K81), or F-type (F21 ) 
may be used to generate an arithmetic node. 

• If area of K21 or thickness of K31 or K41 is input as 0. then a 
surface boundary node is generated. 

• Density (p) must equal 0. 


Standard Conductor - 

. A conduction conductor is generated for all the K21, K31,or K41 elements 
with adjoining edges for all K61 and K81 elements with adjoining 
faces . 

. Surface nodes are considered to have faces only. 

* Negative values of thermal conductivity (K) for temperature 
dependent K. 

• ARRAY tables must be manually input by the user. 


One-Way Conductor - 

• A one-way conductor is generated for all F-type (F21) elements 
having a common edge, connecting the eg of these elements. 

• The conductor value is calculated a Cp A where, M, flowrate 
per unit flow area (input in the K location of the NMAT card) 
Cp, specific heat, and A input cross-sectional area. 

• Negative value of flowrate (fi) for temperature dependent M. 

• ARRAY tables must be manually input by the user. 


Convection Conductor - 

• One standard conductor generated for each C-type (C21 , C31 or C41) 
element input connecting the element defined by the grid points 
input on the C-type element and the element on the NFILM card. 

• NFILM card is required for each different convective sink node 
or convective rate (He). 


. NFILM cards are manually input by the user. 

. Negative value of convective rate (He) for temperature dependent He. 

. ARRAY tables must be manually input by the user. 

. An arithmetic surface node is created on the element to which the 
C-type element is applied except when applied to surface elements 
and 2- grid point elements. 

. Only one surface node is created regardless to the number of C~type 
elements applied to the surface. 


Radiation Conductor - 

. One radiation conductor generated for each R-type (R21, R31, or R41) 
element. The connecting elements are defined by the grid points on 
the R-type element and the element on the NRAD card. 

. NRAD card is reguired for each different radiation sink node or 
radiative rate (eaF). 

. NRAD cards are manually input by the user. 

. Negative value of radiative rate (eaF) for temperature 
dependent eaF. 

. ARRAY tables must be manually input by the user. 

. An arithmetic surface node is created on the element to which the R-Type 
element is applied except when applied to surface elements and 2 grid 
point elements. 

. Only one surface node is created regardless to the number of R-Type 
elements applied to the surface. 
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1.4 MODELING TECHNIQUE. 


To familiarize the user with the method of thermal modeling using 
the finite elements discussed in Section 1.2, Figures 2 thru 4 are 
presented. Figure 2 presents a typical model of a conduction element 
with radiation and convection. An arithmetic surface (K41) is placed 
on the conduction element (K81) so the thermal resistance between the 
conduction element and the surface will be accounted for in the lumped 
parameter model. If desired, the !<41 element could have been excluded. 

A convection element (C41) and radiation element (R41) are applied to 
the surface element. As many C41 and R41 elements required can be 
applied to an element. For this example, the convection and radiation 
elements are connected to a K41 boundary element. These elements 
can be graphically displayed to aid in modal check-out. The lumped 
parameter designation will appear as in Figure 2 with D1 representing 
the K81 element, A2 representing the arithmetic element (K41), and 
B4 representing the boundary element. "D" signifies a diffusion node, 

“A" signifies an arithmetic node, and “B" a boundary node. The 
conductor network. can also be graphically displayed. As in the labeled 
conductors of the example, "K" signifies a standard conduction conductor, 
"C" a convection conductor, and "R" a radiation conductor. 

Figure 3 presents a typical pipe/fluid model. This figure shows 
the finite element representation, the thermal idealization of the 
model, and the graphical display of thermal model. It is important to 
note in fluid flow analysis, that F-type elements must be used because 
they represent thermal nodes which are connected by “one-way" conductors. 
F-type conductors are graphically represented by "F" preceeding the 
conductor number. 

When a transition from 1-D to 2-D or 3-D, or 2-D to 3-D elements 
is required, special modeling is required. Figure 4 gives the method 
required to change from 1-D to 2-D elements. As shown in the figure, 
an arithmetic surface is applied to (1) and (2) and a zero length 
arithmetic surface added to the element (4). These arithmetic surfaces 
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cause conductors K101, K1G2, and K103 to be added to the thermal network, 
but these conductors aren't connected (see graphical display of model ). 

To connect these conductors, nodes A5, A6, and A7 must be converted to 
a single node. That is delete two of these nodes from the MITAS node 
data block and convert these deleted node numbers to the remaining node 
number on all the conductor cards in the MITAS conductor block. 

Resulting in the thermal idealization shown in the figure. These 
changes can be made with the text editor. Other dimensional transitions 
are similarly made. 
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B4 / FREE STREAM TEMPERATURE 

X K41 (BOUNDARY MODE) 

/ ® GRID POINT 

Xo (n) FINITE ELEMENT ID 

C502{ Convection or 
Radiation 

R503 Conductor) 7 C41 (CONVECTION) 

h» or 

/ R41 (RADIATION) 


K41 (ARITHMETIC SURFACE) 



V 
. 6 . 



|K501 (CtyiducMon x K81 (CONDUCTION ELEMENT) 
lni ' Con doctor) 


* NOTE: C41 and R41 will perform function of K41 (arithmetic surface) 
if desired. Above shown for demonstration only. 


Figure 2 . Typical Convection or Radiation Model 
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Figure 3 . Typical Pipe/Fluid Model 





o - DIFFUSION NODE 
X - ARITHMETIC NODE 



USER MANUALLY 
CHANGES A5 -» A6, 

A7 ->A6 USING 
TEXT EDITOR ON 
MITAS CONDUCTOR 
CARDS, AND DELETE 
A5 AND A7 NODE CARDS 


Figure 4. Conduction Transistion From 1-D to 2-D Elements, 
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Section 2 


THERMAL MESH GENERATION -SPAR PROGRAM 

The SPAR program is a finite element structural and thermal analysis 
program which operates on a common data base using mutiple processors. 

Three of these processors are designed to input to the data base* all of the 
necessary physical and material properties of the thermal problem. These 
processors* TAB, TELD and AUS are well suited for the task of finite 
element mesh generation since they utilize free field input and have 
extensive options for automated mesh generation. The TAB arid TELD processors 
have been modified to output joint location and element connection data 
cards of the thermal model, which can be interpreted by the CINGEN program. 

The following sections describe the capability of the SPAR, TAB, TELD 
and AUS processors for mesh generation, illustrate typical deck setups, and 
describe input and output data cards. 

2. 1 CAPABILITY 

Within the TAB and TELD processors are several subprocessors which 
perform specific input and output tasks. The following subprocessors are 
used for thermal model mesh generation. 

Processor TAB (CPC). or TABX (UNIVAC) 

Parameter: 

ONLINE - A parameter used to control printout and punch card output. 

Subprocessors: 

START - Used to define the number of joints or grid points used 

in the finite element model; list of excluded joint motion 
components. 

JLOC - Used to define the joint location (X,Y,Z) with a reference 
coordinate system. This subprocessor allows for automatic 
mesh generation through "loop- limit" logic. 
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Processor TELD (CPC or UNI VAC) 


The TELD processor provides for thermal element definition. Elements 
which may be generated include conduction, one-way flow conductors, con- 
vection and radiation. The subprocessors described below are used to 
select the type of element and to direct card output for the CINGEN 
program. • 

Parameter, 

CINGEN - A RESET parameter used to control punch card output. 

NUTED - A RESET parameter used to select thermal element generation. 

Subprocessors, 

K21 , K31, K41 , K61, K81 - Used to generate 1-, 2-, and 3-dimensional thermal 
conduction finite elements. These subprocessors allow for automatic 
element generation through “Loop-Limit" logic. 

Note: To obtain a one-way conductor (F21), a K21 element must be 

placed in element Group 2 during processing in the TELD 
processor (see Appendix A for a description on placing 
selected elements in a GROUP). 

C21 , C31, C41 - Used to generate 2-, 3-, and 4-node thermal convection 
surfaces. These processors allow for automatic element generation 
through "Loop-Limit" logic. 

R21, R31 , R41 - Used to generate 2-, 3-, and 4-node thermal radiation 
surfaces. These processors allow fcr automatic element generation 
through "Loop-Limit" logic. 

Processor AUS (CPC or UNIVAC) 

The AUS processor is a general matrix and table manipulation package 
which is used to input material and section properties, convection 
coefficients and radiation coefficients when generating a thermal finite 
element model. The procedure for using this processor has not been 
modified and no card output occurs from this processor. SPAR data blocks 
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are generated for use in the TELD processor. Therefore, the AUS processor 
must be used to input properties prior to execution of the TELD processor. 

TABLE - The table processor is used to input the SPAR thermal 

analyzer the conduction, convection and radiation properties, 
conductor thickness, and surface areas. 

2.2 LIMITATIONS 

The SPAR thermal analyzer is still under going development. As a 
result, material properties were not punched out automatically during 
TAB, AUS or TELD execution. These properties must be prepared manually by 
the user for input to the CINGEN program. This was not considered a severe 
limitation since most thermal problems involve only a small number of 
dissimilar materials. 

Temperature dependent conduction, convection and radiation properties 
defined on SPAR input cards COND PROP, etc. are currently not used to 
generate the MITAS (or SINDA) temperature dependent ARRAYS. . SPAR INPUT 
cards are described, however, for completeness; since some users may choose 
to analyze a problem using MITAS and SPAR. 

There are no problem size restrictions when using SPAR for thermal 
model mesh generation. However, when large problems are involved, the 
user must RESET the CORE in the manner prescribed by the SPAR program 
and adequare mass storage must be provided for holding the punch card 
output. 

2.3 INPUT/OUTPUT DESCRIPTION 

This section describes each type of input card available for generating 
grid mesh and thermal finite element networks. The OUTPUT from SPAR is 
punched cards. The output format of each card is described under 
"REMARKS" on each input card description. 
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SPAR DATA CARD DESCRIPTION 


Subprocessor : JLOC 

Description : Use to produce the position coordinates of the finite 

element joints. 

Format and Exampl e: 

JLOC 

K, XAl • XA 2 , XA 3 , XBl, XB 2 , XB 3 , NI, IJUMP, NJ, 

(if UJ is given a second card must appear) 

JJUMP, XCi, XC ? , XC 3 , XD], XD 2 , XD 3 

JLOC 

1, 0.1, 5.2, 3.0 (Figure A) 

66, 1.0, 1.0* 2.0, 10.0, 1.0, 2.0, 10 (FigureB) 

5, 1.0, 1.0, 2.0, 10.0, 1.0, 2.0, 10, 2, 2 (Figure C) 

1 , 1.0, 1.0, 4.0, 10.0, 1.0, 4.0 


Field Contents 

JLOC Subprocessor of TABX (UNI VAC) or SPAR PUNCH TAB 


K 

XAi, XA 2 , XT 3 
XBl, XB 2 , XB 3 
XCl , XC 2 , XC 3 
XDl, XD 2 , XD 3 
NI 

IJUMP 

NJ 

JJUMP 


Joint index 

Joint A coordinate location 

Joint B coordinate location 

Joint C coordinate location 

Joint D coordinate location 

Number of joints defined in the I-direction 

Increment added to the joint index in the I-direction. 

Number of joints defined in the J-direction 

Increment added to the joint index in the J-direction. 
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OPTION 1 


OPTION 2 


Only one joint (k) is 
define at location 
XAi, XA 2 , XA 3 


Point B 

l- ' 

/*K+(NI-1) ( I JUMP) 



K = 1 


• (XAl, = 0-1, 
XA2 = 5.2, 
XA 3 = 3.0) 


/ 

K+2 (I JUMP) 


Point A '/ 

y 

s K 


y.\ 


%K+IJUMP 


Figure A 


Fi gure B 


OPTION 3 


Point D 


TA3X 


Point C 

! 

K + ( J JUMP ) 


® K=(NJ-1 )JJUMP 


Point A 


© 


K 


+ 4 

K + IJUMP + JJUMP 


K + tjUMP 


NI Joints 


K + 2( IJUMP) + 


(NI *NJ joints are defined) 
Figure C 

Remarks : 


© 


NJ Joints 


© 

Point B 
K-(NI-l) IJUMP 


1. There are three possible interpretations of the input card. Each 
is described in Figures A, B, and C. 

2. TABX (Univac) or Sparpunch TAB (CDC) must be used to obtain Grid 
card punch output for use by the CINGEN or GEOMPLT programs. 


3. Punch card output is selected by the ONLINE parameter. Where, 


ONLINE - 0 Minimum printout 

1 Normal printout (DEFAULT) 

2 Maximum printout 

3 Print SPAR output and punch NASTRAN GRID card with 
displacement coordinate system (Field 7) blank. 

4 Print SPAR output and punch NASTRAN GRID card with 
displacement coordinate system (Field 7) equal 

to local grid point coordinate system. 

Example: 

@XQT TABX 
ONLINE = 3 
Start 125 
JLOC 

5,1.0,1.0,10.7 

4. The GRID card punch output format is: 

GRID GID GD XI X2 X3 CD (A4, 4X, 218, 3F8.3, 18) 


Where, 

GID - Grid point number. . 

GD = Local coordinate system used to define GRID location 

XI = Local X position coordinate 

X2 = Local Y position coordinate 

X3 = Local Z position coordinate 

CD » Displacement coordinate system. 

5. See Appendix A for detail discussion on JLOC prodecures. 
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SPAR DATA CARD DESCRIPTION 


S ubprocessor : K21 , C21 or R21 

Description : These subprocessors are used to generate a 2-node, one- 
dimensional conduction, convection or radiation element. 


Format and Example {Typical ): 


K21 


NMAT = nmat 
NSECT * J 

Jl, J2, Netopt, HI, JK, JINC (if Netopt * 1 or 2) 
Jl, 02, Newtop, NI, IINC, NJ, JINC (if Netopt =3) 

K21 

NMAT = 5 
NSET = 2 
NSECT 

9, 12, 1, 4, 2, 100 (Figure a) 

9, 12, 2, 4, 2, 100 (Figure b) 

4, 37, 3, 4, 10, 2, 100 (Figure c) 


Contents 


2-node conduction element subprocessor 
2-node convection element subprocessor 
2-node radiation element subprocessor 

First joint number of first element generated 
Second joint number of first element generated 
Element network option 

Number of elements generated in the I-direction 
Number of elements generated in the J-direction 
Increment added to joints in the I-direction 
Increment added to joints in the J-direction 

Material property table pointer (nmat, nfilrn, nrad) 

Section property table pointer (J) 

NETOPT =1 i NETOPT * 2 NETOPT * 3 
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Field 

K21 

C21 

R21 

Jl 

J2 

Netopt 

NI 

NJ 

IINC 

JINC 

NMAT 

NSECT 


Remarks: 


1. The order In which the elements are defined is indicated by the number 
enclosed in parentheses. The index number (EID) identifying elements 
within each group are determined by the order in which the elements 
are defined. 

2. Conduction properties are selected by the table pointer NMAT = nmat, where 
nmat = table number {default = 1) of the data block "CORD PROP NMAT 9". 

3. Section properties (Area) are selected by the table pointer nsect = J, 
where 0 = table number (default = 1) of the data block "K Area nsect 9" 
for conduction elements * : 'C CIRC nsect 2" for convection element; 

and "R CIRC nsect 2" for radiation elements. 

4. Convection properties (h) are selected by the table pointer NMAT = nfilm 
where nfilm “ table number (default = 1) of the data block "CONV PROP nmat 
2 ". 

5. Radiation properties (eoF) are selected by the table pointer NMAT = nrad, 
where nrad = table number (default = 1) of the data block" RADI PROP nmat 
2 ". 

6. One card is punched for each element if the parameter CINGEN = 1 is 
RESET. Punch card format is, 

K21 EID nmat J1 J2 Area (A3, 5X, 418, 1PE16.7) 

C21 EID nfilm J1 J2 Acirf, (A3, 5X, 418, 1PE16.7) 

R21 EID nrad J1 J2 Acirf. (A3, 5X, 418, 1PE16.7) 

7. A single element may be generated by specifying only Jl, J2. 

8. Surface elements are defined for CINGEN arithmetic nodes if: 

Area =0.0 for K21 elements 


9. Arithmetic nodes are defined at the centroid of C21 and R21 elements. 
If more than one element type exists at the same spatial location, 
CINGEN will generate only one MITAS arithmetic node. Nodes are not 
created for 2-grid point elements. 

10. To obtain a one-way conductor (F21), a K21 element must be placed in 
element GROUP 2, e.g. , 

@XQT TELD 

RESET CINGEN = 1 , NUTED = 1 
K21 

GROUP T Pipe Conduction ) 

46 47 1 31 ) Punch K21 elements 

50 51 1 4 1 ) 

GROUP 2' One-way Flow Conduction ) 

55 56 1 21 ) Punch F21 elements 


SPAR DATA CARD DESCRIPTION 


Subprocessor: K31, C31, R31 

K41 , C41 , R41 

Description : These subprocessors are used to generate 3- or 4-node two- 

dimensional conduction, convection or radiation elements. 

Format and Example (Typical): 

K31 

NMAT = nmat 
NSECT = J 

Jl, J2, J3, Netopt, NI, NJ (if Netopt = 1 

Jl, 02, J3, Netopt, NI, NJ, JINC ( if Netopt = 2 or 3) 

K41 

NMAT = nmat 

Jl, J2, J3, J4, Netopt, NI, NJ, NK, KINC (if Netopt = 1) 

Jl, J2, J3, J4, Netopt, NI, NJ (if Netopt = 2) 

K31 

NMAT = 5 
NSECT =3 

2, 3, 103, 1, 3, 2 
2, 5, 7, 2, 6, 30 
2, 5, 7, 3, 6, 2, 30 
K41 

2, 3, 23, 22, 1 , 2, 3, 2, 200 (Figure d) 

1, 11 12, 2, 2, 6, 3 (Figure e) 


Contents 

3- or 4-node conduction element subprocessor 
3- or 4-node convection element subprocessor 
3- or 4-node radiation element subprocessor 
first joint number of first element generated 
second joint number of first element generated 
third joint number of first element generated 
fourth joint number of first element generated 
element network option 

number of elements generated in the I-di recti on 
number of elements generated in the J-di recti on 
number of elements generated in the K-direction 
increment added to joints in the K-direction 
material property table pointer (nmat, nfilm, nrad) 
section property table pointer (J) 


Field 

K31 ,K41 
C31 ,C41 
R31 ,R41 
Jl 
J2 
J3 
J4 

Netopt 

NI 

NJ 

NK 

KINC 

NMAT 

NSECT 


(Figure a) 
(Figure b) 
(Figure c) 
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I EiL-LF 


3 NODE, NETOPT * 1 
202 203 204 205 


102 


/ 103 

(To V 

/ 104 

(12)/ 
/)1) 
/ 105 

(21/ 

/fn 

(4)/ 

>(3) 

l>(5) 


Figure a 


3 NODE, NETOPT = 2 


3 NODE, NETOPT = 3 



Figure c 


4 NODE, NETOPT = 1 



Figure d 

Remarks : 


4 NODE, NETOPT = 2 



32 

31 


1. The order in which the elements are defined is indicated in the above 
figures by the number enclosed in parentheses. The index number (EID) 
identifying elements within each group are determined by the order in 
which the elements are defined. 


2. Conduction properties are selected by the table pointer NMAT= nmat, where 
nmat = table number (default = 1) of the data block COND PROP nmat 9". 


3. Section properties (chicxness) are selected by the table pointer NSECT=J, 
where J = table number (default = 1) of the data block "K THIC nsect 2" 
for conduction elements. 

4. Convection properties (h) are selected by the table pointer NMAT = nfilm 
where nfilm - table number (default = 1) of the data block, "CONV PROP 
nfilm 2". 

5. Radiation properties (eoF) are selected by the table pointer NMAT «* nrad 
where nrad = table number (default =1) of the data block "RADI PROP 
nrad 2". 
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SPAR DATA CARD DESCRIPTION 


Remarks : 

6. One card is punched for each element if the parameter CINGEN = 1 is RESET. 
Punch card format is, 

K31 EID HID J1 J2 J3 thick (A3, 5X, 518, 1PE16.7) 

C31 EID NFILM J1 «J2 J3 (A3, 5X, 518) 

R31 EID NRAD J1 J2 J3 (A3, 5X, 518) 

K41 EID MID ul J2 J3 J4 thick (A3, 5X, 618, 1PE16.7) 

C41 EID NFILM J1 J2 J3 J4 (A3, 5X, 618) 

R41 EID NRAD J1 J2 J3 J4 (A3, 5X, 618) 

7. A single element may be generated by specifying only J1 , «J2, J3, or 
dl, J2, J3, J4. 

8. Surface elements are defined for CINGEN arithmetic nodes if: 

Thick =0.0, for K31 and K41 elements. 

9. Arithmetic nodes are defined at the centroid of C31 , C41 , R31 , and R41 
elements. If more than one element type exists at the same spatial location, 
CINGEN will generate only one MITAS arithmetic node. 
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TELD 


SPAR DATA CARD DESCRIPTION 


SUBPPQCESSOR: K61 

Description : Subprocessor used to generate a 6-node pentahedron conduction 

element. 

F ormat and Example : 

K61 (A3) 

NMAT = nmat 

JA a JB, NJ, NK, JINC, KINC, ICLOSE, NETOPT (free field, use space or 

, comma as delimiter) 

K61 

NMAT * 5 

1, 2, 5, 1 , 1 , 200, 1 ,1 


Field 


Contents 


K61 


JA 

JB 

NJ 

NK 

JINC 
KINC 
I CLOSE 

NETOPT 

NMAT 


Element generation processor for 6-node pentahedron conduction 
element 

First node of a fan of elements (pivotal node) 

Second node (first node on outer edge of fan) 

Number of elements in the fan (default =1) 

Number of element layers (default =1) 

Increment added to JB resulting in index for JC = JB + JINC 
Increment added to JA, JB, JC resulting in KA, KB, KC 

0 = Open Fan (default) 

1 = Closed Fan 

Element Network option, 0 = Single element (default) 

1 = Fan network 

Table pointer for conduction properties (default =1) 

See COND PROP input card description 



CLOSED PENTAHEDRON FAN 
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REMARKS : 


1. Conduction properties are selected by the table pointer, 

NMAT = nmat, where nmat = table number (default = 1) of the data block 
"COND PROP nmat 9." 

2. Element identification numbers (EID) are assigned by the SPAR Program. 

3. One card is punched for each element if the parameter CINGEN is RESET 
to 1 at the start of the TELD processor e.g. , 

0XQT TELD 

RESET CINGEN = 1, NUTED = 1 
K61 

etc 

4. Punch card format is: 

Card 1: K61 EID MID JA JB JC KA KB KC (A3, 5X, 818) 

5. When NET0PT = 0, input 6 joints -for a single element. 


•3 
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i CLU 


SPAR DATA CARD DESCRIPTION 



Subprocessor : K81 

Description : Subprocessor used to generated an 8-node hexahedron conduction 

element. 


Format and Example: 


K81 

JO, NI, NJ, NK, I INC, JINC, KINC, ICLOSE, NETOPT 


K81 

NMAT=2 

1, 1, 2, 3, I, 100, 10, 0, 1 


Field Contents 

K81 Element generation subprocessor for 8-node hexahedron 

conduction element. 


JO Starting joint number 

NI Number of elements in the I-directian 

NJ Number of elements in the J-direction 

NK Number of elements in the K-direction 

IINC Increment added to joints in I-direction 

JINC Increment added to joints in J-direction 

KINC Increment added to joints in K-direction 

ICLOSE 0 = OPEN string of elements in I-direction 

1 = CLOSED (Ring) string of elements in I-direction 
NMAT Table pointer for conduction properties (default = 1) 

NETOPT Element network option 0 = single element (default) 

1 = network in 3 directions 





TELD 


REMARKS : 

1. Conduction properties are selected by the table pointer, NMAT = nmat 
where nmat = table number (default =1) of the data block “COND PROP 
nmat 9." 

2. Element identification numbers (EID) are assigned by the SPAR Program 
in the order they are generated, 

3. Two cards are punched for each element if the parameter CINGEN is RESET 
to 1 at the start of the TELD processor. 

4. Punch card format is: 

Card 1: K81 EID NMAT JA JB JC JD KA KB (A3, 5X, 818) 

Card 2: KC KD (8X, 218) 

5. When NETOPT = 0, input 8 joints for a single element. 
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SPAR DATA CARD DESCRIPTION 


Table: K AREA 

or 

K THIC 


Descriptio n: Defines tables of conduction area or thickness 

Format and Example : 

TABLE (NI = 1, NJ « j) •, K AREA 1 ; J = 1, j 
Ai » A2> * * A j 

TABLE (NI = 1 , NO = 3) i K AREA 1 ; 0 = 1 , 3 

0.16,2.45,1.73 


OR 

TABLE (NI * 1, Nd = j) 5 K THIC 1 ; J * 1, j 
t] , tg ...... tj 

TABLE (NI = 1, NO = 1) ; K THIC 1 ; J = 1 
0.025 


Field 


Contents 


TABLE 

$ 

A 

t 


name of the AUS subprocessor used to create a new table 
in the SPAR data complex. 

number of areas or thickness input in table. 

area 

thickness 


Remarks : 

1. The thickness is selected by the table pointer nsect = J, during 
K31 or K41 conduction element generation in processor TELD. 

2. The area is selected by the table pointer nsect = J, during K21 
Conduction element generation in processor TELD. J points to a data 
item in "K AREA .1" Table. 

3. No punch card output is created for CINGEN. 

4. The 5/8 punch (: on UNIVAC) or 12/8/7 punch (; on CDC) terminates one 
record and initiates a new record on the same card. 
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AUS 


SPAR DATA CARD DESCRIPTION 


Table: C CIRC 

R CIRC 

Description : Defines tables for convection and radiation surface areas for 1 

dimensional elements. 

Format and Example: 


TABLE (HI » 1. NJ - j) 

c 1 # c 2 ,... Cj 

; c 

CIRC 

1 

; J = 1 

TABLE (NI = 1 , NJ “ 3) 
0.16, 2.45, 1.73 

; c 

CIRC 

1 

; J = 1 

OR 





TABLE (NI = 1, NJ = j) 
C-j, C2» ...Cj 

; R 

CIRC 

1 

5 J - 1 

TABLE (NI — 1 * NJ = 1) 

; R 

CIRC 

1 

; J “ 1 


0.025 

Field 

TABLE 


Contents 

name of the AUS subprocessor used to create a new table 
in the SPAR data complex 


J 

C 


number of circumferences 
circumference 


input in table 


Remarks: 

1. The circumference is selected by the table pointer ,nsect = J, during 
C21 or R21 convection or radiation element generation in processor 
TELD. 

2. No punch card output is created for CINGEN. 

3. The 5/8 punch (: on UNIVAC) or 12/8/7 punch (; on CDC) terminates one 
record and initiates a new record on the same card. 
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AUS 


SPAR DATA CARD DESCRIPTION 


Table: COND PROP 

Description : Defines the temperature dependent property tables. 

Format and Example : 

Table (NI = 1, NJ = j) ; COND PROP nmat ; I = 1, 2, .. NIS d * 1, j 
T, p, c, i s K?2* K 33* K 12* K 23* ! Sl 

Table (NI = 9, NJ » 1 ) i COND PROP 1 ; I* 1 , 2, 3, 4 i J=1 

0 ., 1 ., 1 ., 20 

Field 

TABLE 

i 

j 

nmat 
T 

P 


Remarks : 

1. Conduction properties are selected by the table pointer, NMAT = nmat 
during conduction element (K21, K31 , etc) generation in processor TELD. 

2. If conduction properties are temperature dependent, set NJ = j where 

0 equal the number of temperatures used to define the material properties, 
"j" data lines follow the TABLE card, each data line defines the properties 
corresponding to Temperature T. 

3. No punch card output is created for CINGEN- 

4. All required K-. must be input. 

'J 


Contents 

name of the AUS subprocessor used to 
create a new table in the SPAR data complex 

number of constants input 

number of data lines corresponding to temper 

ature T. 

material ID (integer > 0) 

Temperature 

mass density (at temperature T) 
specific heat (at temperature T) 

conductivity (at temperature T)defined with 
respect to the element coordinates. Default 
= 0 . 
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SPAR DATA CARD DESCRIPTION 


Table : CONV PROP 

DESCRIPTION : Defines the temperature dependent convective film 

coefficient. 

Format and Example : 

TABLE ( NX = 2 NJ = j) ; CONV PROP nfilm ; J = 1, j 
T, h 

TABLE (NI = 2, NJ - 2) ; CONV PROP nfilm; J = 1, 2 

100., 1. + 5 

200. , 1 .5 + 5 

Contents 


name of the AUS subprocessor used to create a new table in 
the SPAR data complex. 

number of data lines following the TABLE card. 

film coefficient property ID (integer >0) 

T Temperature 

h film coefficient 

Remarks : 

1. Convection properties are selected by the table pointer, NMAT = nfilm, 
during convection (C21, C31, C41) element generation in processor TELD. 

2. If convection properties are temperature dependent, set NJ = j, where 
j equal the number of temperatures used. 

3. No punch card output is created for CINGEN. 


Field 

TABLE 

j 

nfilm 


AUS 


SPAR DATA CARD DESCRIPTION 


Table: RADI PROP 


Description : Defines the temperature dependent radiation coefficient. 

Format and Example : 

TABLE (NI = 2 NJ = j) ; RADI PROP nrad 5 J = 1, j 
T, £ 


TABLE (NI = 2, NJ = 2) ; RADI PROP nrad ; J = 1, 2 

100., 0.5 

200 ., 0.6 


Field 

TABLE 

0 

nrad 

T 


Contents 

name of the AUS subprocessor used to create a new table 
in the SPAR data complex. 

number of data lines following the TABLE card. 

radiation coefficient ID (integer >0) 

Temperature 

emissivity 


Remarks 


1. Radiation properties are selected by the table pointer, NMAT = nrad, 
during radiation (R21, R31, R41) element generation in processor TELD. 

2 . If radiation properties are temperature dependent, set NJ = j, where 
j equals the number of temperatures used. 

3. No punch card Output is created for CINGEN. 
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Section 3 

FINITE ELEMENT TO LUMPED PARAMETER CONVERSION - CINGEN PROGRAM 

CINGEN is a computer program developed by Sperry Support Services 
to aid the analyst in the development of thermal models. The program 
utilizes finite element input data, which can be generated from mesh 
generators, such as the SPAR TAB and ELD processors. This data is 
used to create a lumped parameter themal network formatted for thermal 
analyzer programs, such as CINDA, SINDA, and MITAS. CINGEN also creates 
data files for the production of temperature graphs and thermal network 
plots by the GEOMPLT program. 

3,1 CAPABILITY AND METHOD 

The CINGEN computer program has the capability of accepting SPAR or 
NASTRAN finite element data and generating l-» 2-, and 3-dimensional 
lumped parameter thermal networks. Permissable element types include: 
conduction, one-way flow conductors, convection, and radiation. The 
conduction elements can have 2, 3,4, 6, or 8 grid points, the one-way 
flow element can have 2 grid points and radiation and convection elements 
can have 2,3, and 4 grid points. In each conduction or one-way flow, 
element, a thermal node is assigned to the element. These nodes can be either 
diffusion, arithmetic, or boundary nodes as determined by the density value 
selected for the element. Arithmetic and boundary type surface nodes are 
also available. Figure 5 presents a functional flow diagram of CINGEN. 

In determining a diffusion node capacitance (pC v V) a positive, real value 
must be input for the density (p). The specific heat (C v ) can be a positive, 
real value or a negative, real value for temperature dependent property. 

The volume (V) is calculated by the CINGEN program for each finite 
element. The volume of six (6) and eight (8) grid point elements is 
calculated by dividing the element into tetrahedrons and employing vector 
calculus to yield the exact volume. The volume of three (3) and four (4) 
grid point elements is calculated as the product of the element surface 
area (A) and the element thickness (t). The area is computed from the 
coordinate location of the element grid points, and the thickness is an' 
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Figure 5, CINGEN Functional Flow Diagram 
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Figure 5 


CREATE FILES 
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Figure 5 (continued) 
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input parameter (see CINGEN Data Card Section). The volume of a 2- 
grid point element is calculated as a product of the length of the 
element and the cross sectional area. The length is determined by the 
coordinate location of the grid points and the cross sectional area is 
an input parameter (see CINGEN Data Card Section). If there is variable 
specific heat (temperature dependent), the diffusion node would be output 
by CINGEN as a single interpolated node (MITAS type CGS). Variable 
specific heat is specified by using a negative specific heat on the 
element material card. The absolute value of the negative specific 
heat is used as the specific heat array number on the MITAS CGS card. 

(User must supply MITAS ARRAY DATA BLOCK.) 

The center of gravity of the element is computed. Thermal resistance 
from the center of gravity to each face or edge is determined by the equation. 



where 

K = thermal conductivity 

1^= distance from eg of element to the eg of the surface or edge 

Vol = volume of element. 

This equation produces accurate values of thermal resistance for 
elements with constant cross-seci tonal area perpendicular to the direction 
of heat flow. Finite element modeling techniques require that element 
(triangular and quadrilateral plates) aspect ratios less than 3 to 1 be 
used for structural acceptable accuracy. This same rule should be applied 
to thermal modeling to assure modeling accuracy. 

Temperature dependent thermal conductivity is defined similar to the 
variable specific heat. That is the conductivity is input as a negative 
value. Each edge (ID and 2D element) or face (3D element) is then assigned 
a coded number. After all of the conduction and one-way elements have been 
processed, the coded numbers are compared to determine common edges or 
faces. When two elements have a common edge or face, they are connected 
by a conductor equal to the reciprocal of the sum of the resistance to each 
face. 
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After all the conduction and one-way elements are connected by 
conductors the convection elements are processed. Convection can be applied j 

to any type of element. The grid points on the convection element card are j 

assigned a coded number then compared with the coded numbers for the con- 
duction and one-way elements. When common edge or face is found a conductor 
is connected to that element given on the NFILM card. The convection rate 
(he) is defined on the NFILM card (can be temperature varying) and the area ] 

is detrained from the grid points on the convection card or read from the __ j 

element card for 2 grid point nodes. Radiation conductors are processed after j 

convection conductors. Radiation is applied the same way as convection ex- j 

^ cept NRAD cards are used to define the radiative rate (eaF). CINGEN multiplies 1 

this value by the area. An arithmetic node is generated on the surface. i 

to which convection and radiation are applied. Only one node is generated 
for each surface. 

3.2 LIMITATIONS 

CINGEN has built-in program limitations which were used to reduce core 
requirements. These limitations are the number of material property cards and 
number of nodes. Explanation of these limitations and how to increase the 
number of material property cards and nodes is given below. 

1. NMAT cards «u«f®*itly limited to 50 which should be adequate for 

most therma': If more NMAT cards are required the following 

dimensional variants in common block /PROP/ of CINGEN (main program and 
subroutine SPAR) must have their dimensions increased to the desired 
limit. These variables are MAT7, RHO, CONDUC, CP, CV, CFLAG, and CVFLAG. 

2. CINGEN is currently limited to a maximum of 500 nodes. This maximum 
can be changed by changing the following cards in the main routine of 
CINGEN 

COMMON BUCKET (19500) 

LOAN = 19500 

The number of nodes allowed is determined by the core size assigned 
(19500) to COMMON BUCKET. For 1000 nodes these cards would be 

COMMON BUCKET (39000) e.g. CORE _ NODES 

19500 ~ 500 

LOAN = 39000 

3. The element ID (EID) and CINGEN node number must be the same if a 
node and conductor network are to be graphically displayed by GEOMPLT. 

That is, NODADD must be input as 0 (See page 3-13). 
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3.3 INPUT/OUTPUT DESCRIPTION 

CINGEN is capable of processing both SPAR and NASTRAN elements. 

The acceptable SPAR elements are conduction (K21 , K31 » K41 » K61 and 
K81), one-way flow element (F21), convection elements (C21, C31* and 
C41), and radiation elements (R21* R31» and R41). The acceptable 
property cards are NMAT.» NFILM, and NRAD. NASTRAN type elements are 
CBAR, CTRI , CQUAD, CWEDGE, and CHEX. The program keys only on the 
first two characters of these element names. CINGEN reads data from 
files as assigned by the variable IMP in the MAIN program. This value is 
currently set to 4 for reading file 4* change INP = 5 for card reader. 
Figure 6 shows the basic data set-up for CINGEN. The data is broken 
into two blocks. BLOCK 1 contains the NMAT 9 NFILM, and NRAD cards. 

The cards can be mixed and in any order. BLOCK 1 must be preceeded 
by '1st Input Record' and ended with an END card. These two cards 
are required by CINGEN. BLOCK 2 contains the TITLE card, GRID cards s 
and ELEMENT cards. These cards can be mixed and in any order. The 
TITLE card is not required. This data block is ended with a 
ENDDATA card. 

CINGEN outputs a complete MITAS data deck except for the ARRAY and 
CONSTANTS data cards. The user must add these cards according to the array 
numbers assigned on the NMAT, NFILM, and NRAD cards and the constants re- 
quired by his particular problem. Figure 7 is an example CINGEN output. 
CINGEN also creates a data file 29 (TAPE 29) which is used by GEOMPLT to 
graphically display the model. Figure 8 defines the variables output on 
TAPE 29 and their formats. 

Data is input in any consistent set of units and the following pages 
describe the format of these data cards. 



OQOOOl 

0 00 

BCD 

6TITLE 

EHAEJPLE 

NUMBER 2 

000002 

000 

BCD 

6 



000003 

000 

END 




000004 

0 00 

BCD 

3N0"E 1 

DATA 


O00C05 

0 00 

CGS 

It 

75*00 * 

All* .6312+01 

000006 

000 

CGS 

2 * 

75*00* 

All* *6312+01 

000007 

000 

CGS 

3* 

75*00* 

All* *6312+01 

000006 

000 

CGS 

4* 

75*00* 

All* .6312+01 

000009 

000 

CGS 

5» 

75.00* 

All* .6312+01 

000010 

0 00 

CGS 

6? 

75.00* 

All* .6312+01 

000011 

0 00 

CGS 

7, 

75.00* 

All* *6312+01 

000012 

000 

CGS 

8* 

75.00* 

All* .6312+01 

000013 

000 

CGS 

9 * 

75.00* 

All* *3156+01 

000014 

ODD 

CGS 

10* 

75*00. 

All* .3156+01 

000015 

0 00 

CGS 

11* 

75*00* 

All* .3156+01 

000016 

0 00 

CGS 

12* 

75.00* 

All* .3156+01 

000017 

000 

CGS 

13* 

75.00* 

All* .6312+01 

000018 

000 

CGS 

14* 

75.00* 

All* .6312+01 

000019 

000 

CGS 

15. 

75.00* 

All* .6312+01 

000020 

000 

CGS 

16* 

75.00* 

All* .6312+01 

000021 

□ 00 

CGS 

17* 

75.00* 

All* .6312+01 

000022 

0 00 

CGS 

18* 

75.00* 

All* .6312+01 

00002 3 

0 00 


47 * 

75.00* 

*0000 

000024 

000 


46* 

75.00* 

.1025+00 

000025 

0 00 


49* 

75.00* 

.1025+00 

000026 

000 


56* 

75.00* 

.1300-01 

000027 

000 


57* 

75.00* 

.1380-01 

000028 

000 


19* 

75.00* 

-.1000+01 

000029 

000 


20* 

75.00* 

-.1000+01 

000030 

000 


21* 

75.00* 

-.1000+01 

000031 

000 


22* 

75.00* 

-.1000+01 

000032 

000 


23* 

75.00* 

-.1000+01 

000033 

oao 


32* 

75.00* 

-.1000+01 

000034 

000 


33? 

75.00* 

-.1000+01 

000035 

000 


34 * 

75.00* 

-.1000*01 

000036 

000 


35* 

75.00* 

-.1000+01 

000037 

0 00 


36* 

75.00* 

-.1000+01 

000038 

000 


-50* 

75.00* 

.0000 

00QQ39 

000 


-51* 

75.00* 

.0000 

000040 

000 


-52* 

75.00* 

.0000 

000041 

0 00 


-53* 

75.00* 

.0000 

000042 

0 00 


-54* 

75.00* 

*0000 

000043 

0 00 


-55* 

75.00* 

•0000 

000044 

000 

END 




000045 

000 

BCD 

3 CONDUCTOR DATA 


000046 

000 

CGS 

1* 

19* 

5* A13* .2000+02 

000047 

000 

CGS 

2* 

20* 

6* A 13* .2000 + 02 

000048 

000 

CGS 

3* 

21* 

7* A13 * .2000+02 

000049 

ooo 

CGS 

4* 

22* 

B* A13 * .2000+02 

000050 

000 

CGS 

5. 

23* 

18+ A13* .5625+01 

000051 

000 

CGS 

6 * 

32* 

11+ A13+ .1000+02 

000052 

000 

CGS 

7* 

33* 

12* A13+ .1000+02 

000053 

ooo 

CGS 

8* 

34* 

14* A13+ .2000*02 

000054 

000 

CGS 

9* 

35* 

16* A 13 * .2000*02 

000055 

ooo 

CGS 

10* 

36* 

18* A 13 + .2000+02 


DIFFUSION NODES 


ARITHMETIC NODES 


BOUNDARY NODES 



Figure 7 


Thermal Model Created By CINGEN For MITAS 
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000056 

000 

CGS 

U* 

1* 

2* 

A 13, 

.2500*01 

000057 

000 

CGS 

22* 

U 

3* 

A 13 , 

.2500*01 

000058 

000 

CSS 

13 r 

1* 

5* 

A 13 * 

.1000*02 

000059 

0 00 

CGS 

14, 

2* 

4* 

A 13* 

.2500*01 

000060 

000 

CSS 

25, 

2* 

6* 

A13 * 

.1000*02 

000061 

000 

CGS 

26* 

3* 

4, 

A13, 

.2500*01 

000062 

000 

CSS 

27* 

3* 

7* 

A 13, 

.1000*02 

000063 

000 

CGS 

18* 

3* 

9* 

A 13 , 

.1800*01 

00006# 

oco 

CGS 

19, 

4* 

B* 

A 13* 

•1000*02 

000065 

0 00 

CSS 

20* 

4* 

20* 

A 13, 

.2647*02 

000066 

000 

CGS 

21* 

5* 

6, 

A13* 

.2500*01 

000067 

000 

CGS 

22* 

5* 

7* 

A 13, 

.2500*01 

000060 

000 

CGS 

23* 

6* 

8* 

AI3 * 

•2500*01 

000069 

000 

CGS 

24* 

7* 

8, 

A13, 

.2500*01 

000070 

000 

CGS 

25* 

7. 

11* 

A 13, 

.1800*01 

000071 

000 

CSS 

26* 

Sr 

12* 

A13, 

•2647*01 

000072 

000 

CGS 

27* 

9* 

10* 

A13, 

•1406*01 

000073 

000 

CSS 

28* 

9* 

Hr 

A 13 * 

.5000*01 

00007# 

000 

CGS 

29* 

9* 

13* 

A13 , 

.2812*01 

000075 

000 

CGS 

30* 

10* 

12* 

A 13, 

.5000*01 

000076 

000 

CGS 

31* 

20* 
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Figure 7. (continued) 
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Figure 8. Typical Data File ?9 as Created by CJ.MGEN 
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Card: 1 st INPUT RECORD 


Description : This card increments all element and conductor numbers, 
specifies the type of coordinate system used, and defines initial reference 
temperature of all thermal nodes. 


Format and Example : 


NODADD, CONADD, ICS, T. ... ,, ETYPE 

initial 

0 500 R 75. 


(2I8,7X,A1,F8.3) 


Field 

NODADD 

CONADD 

ICS 

^initial 


Contents 

Integer that will be added to element number to form 
the node number, NODE NO. ® EID+NODADD. 

The conductors are numbered sequently from beginning 
with CONADD. 

^C retangular coordinate system 
=C cylinderical coordinate system 

Initial reference temperature used on all MITAS node cards. 


ETYPE = NAST, NASTRAN type elements 

f NAST, SPAR type input 


4+CONADD 


EID4+N0DADD 


3+CONADD 


EID1+NODADD 


1+C0NADD 
jt iAAA 



2+C0NADC 


EID3+N0DADD 


EID2+N0DADD 


Where, the EIDi are element 
numbers assigned during element 
generation in the SPAR program. 
Then, 

NODE NO. = EID + NODADD 


Remarks : 

1. The EID's shown in the above figure are the element numbers assigned 
during finite element generation in the SPAR program. 

2. The node number assigned by CINGEN will be 

NODE NO. = EID + NODADD 

3. The conductors shown in the above figure are numbered by the 
CINGEN program. The numbering sequence begins with the first 
conductor generated by CINGEN and is numbered in the following way 

1st conductor N0.= 1 + CONADD 

2nd conductor N0.= 2 + CONADD 
* 

Nth conductor N0.= N + CONADD 
3-13 


4. NODADD must equal 0 for plotting (See page 3 - 6 ) 

5. This card must be input. 
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Field Contents 

NMAT NMAT to signal material properties 

MID material ID 

RHO density 

= 0 elements with this MID will be boundary nodes 

<0 elements with this MID will be arithmetic 
surface nodes 

>0 elements with this MID will be diffusion nodes 
and RHO will be their density 

CV specific heat at constant volume 

<0 temperature variable specific heat* absolute 
value of CV used as MITAS array number 

>0 specific heat of element 

COND thermal conductivity 

<0 temperature variable thermal conductivity 

absolute value of COND used as MITAS array number 

>0 thermal conductivity of element 

CP specific heat at constant pressure (used only in one-way flow conductors) . 

If CP is to be temperature dependent put the array number containing the = 
product of the flow rate per unit surface area and CP (ftCP) -j n the COND ■ 
location and 1.0 in the CP location. 

Remarks : 

1. This data card must be manually prepared by the user and inserted on 
file 4 along with the SPAR generated finite element data. 

2. This card should be placed behind the "1 st INPUT RECORD" card and ahead of 
the finite element connection cards on file 4. 
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CINGEN OATA CARD DESCRIPTION 


Card : NFILM 

Descri ption : Inputs convective heat transfer coefficient 
Format and Example : 

NFILM FID ELEM H c {4A6,free field) 

NFILM 500 99 0.1 

Contents 

NFILM to signal convection properties 
convection property ID 

element number to which the H c will be connected 
heat transfer coefficient 

<0 temperature dependent coefficients, 
absolute value of H used as array 
number c 

>0 heat transfer coefficient per unit area 

Remarks : 

1. Conductor value is calculated as node area times H c , 

2. Convection property ID number must be unique with respect to all other 
convection property ID numbers. 

3. The field length for H c is 15 columns. 


Field 

NFILM 

FID 

ELEM 



CINGEN DATA CARD DESCRIPTION 


Card : NRAD 

Description : Inputs Radiation heat transfer coefficient 

Format and Example : 

NRAD RID ELEM SCRIPT (3A8, AT 6, free field) 

NRAD 6 16 Q.2E-10 

Field Contents 

NRAD NRAD to signal radiation properties 

RID radiation property ID 

ELEM element number to which SCRIPT is connected 

SCRIPT eoF 

<0 temperature dependent absolute value of 
SCRIPT used as array number. 

>0 radiative heat rate per unit area 

Remarks : 

1. Conductor value is calculated as node area times SCRIPT. 

2. Radiation property ID numbers must be unique with respect to all other 
radiation property ID numbers. 

3. The field length for SCRIPT is 16 columns. 
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C INGEN DATA CARD DESCRIPTION 


Card : PBAR 

Description : Property card for NASTRAN 2 grid point element (CBAR) 

Format and Example : 

PBAR PID MID A (4A8,free field) 

PBAR 6 9 .47 

Field Contents 

PBAR PBAR to signal material property card 

PID property ID 

MID material ID 

A cross-sectional area 

Remarks : 

1. Property ID numbers must be unique with respect to all other PBAR 
property numbers. 
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CINGEN DATA CARD DESCRIPTION 


Card : PTRIA2 and PQUAD2 

Description: Property card for NASTRAN 3 and 4 grid point elements (CRTIA2 and 

CQUAD2) 


Format and Example: 

PTRIA2 

PQUAD2 

PID HID T BLANK PID 

PID MID T BLANK PID 

PTRIA2 

100 26 1.8 59 

Field 

Contents 

PTRIA2 

PTRIA2 to signal material 

PQUAD2 

PQUAD2 to signal material 

PID 

property ID 

MID 

material ID 

T 

thickness of element 

BLANK 

blank field 


HID T (8A8» free field) 

HID T (8A8, free field) 

96 .74 

property card 
property card 


Remarks : 

1. One or two element properties may be defined on a single card. 

2. PTRIA2 property ID numbers must be unique with respect to all other 
PTRIA2 property numbers. 

3. PQUAD2 property ID numbers must be unique with respect to all other 
PQUAD2 property numbers. 
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CINGEN DATA CARD DESCRIPTION 


Card : END 

Description: This card ends the first data block input. 
Format and Example : 

VALUE 

END 

Field Contents 

VALUE END to denote end of first data block. 

Remarks : 

1. This data card must be input. 


(A8,free field) 
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CINGEN DATA CARD DESCRIPTION 



Card : TITLE 

Description: Title of the model 



TITLE NAME 

TITLE Example problem 2 


(A8, 8A8) 


Field Contents 

TITLE TITLF - signals a title card 


NAME up to 64 characters, model description 


Remarks : 

1. This card supplies title information for the MITAS data deck. 

2. Format of the MITAS title card punched by CINGEN is 
BCD 6( first 36 characters on the title card) 

BCD 6( last 36 characters on the title card) 

3. Every character in the first 72 columns is printed. 

4. This data card may be omitted. 




CINGEN DATA CARD DESCRIPTION 


Card : GRID 

Description : Defines the coordinate locations of the grid point 
F ormat and Example : 

GRID ID 8 X Y Z (6A8 .free field) 

GRID 102 5. 3. 9. 

Field Contents 

GRID GRID to signal grid point card 

ID grid point ID 

B blank field 

X X coordinate of grid point 

Y Y coordinate of grid point 

Z Z coordinate of grid point 
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C INGEN DATA CARD DESCRIPTION 


Card: K21, K31 , K41, K61 or K81 

Description : Conduction elements with 2, 3, 4, 6 and 

Format and Example : 


K21 

EID 

MID 

01 

0,2 ' 

AREAC 


K31 

EID 

MID 

01 

02 

03 

THICK 

K41 

EID 

MID 

01 

02 

03 

04 

THICK 

K61 

EID 

MID 

01 

02 

03 

0.4 

05 06 

K81 

EID 

MID 

01 

02 

03 

04 

05 06 


07 

08 






K31 

42 

13 

8 

9 

10 




Field 

K21 

K31 

K41 

K61 

K81 

EID 

MID 

01, 02, 03, 04, 
05, 06, 07, 08 

AREAC 

THICK 


Contents 

2- node conduction element 

3- node conduction element 

4- node conduction element 
6-node conduction element 
8-node conduction element 
element ID 

material ID 

grid points defining the element 

cross-sectional area 
thickness of the element 
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grid points 


(9A8,free field) 


(9A8,8X,2A8,free field) 


Remarks: 


1. The grid points of K31 and K41 elements must be input in a consecutive 
order around the element. 

2. Grid point numbering of K61 and K81 elements must be input in a 
consecutive order around any face but the opposite face has to be 
input in a similar way starting on the same edge as the first 
face as shown in the following examples: 




3. K21 s K31, and K41 become surface node elements if AREA = 0.0 
for K21 and THICK =0.0 for the K31 and K41 elements. Use 
K31 or K41 elements to obtain surface nodes on faces of K61 
and K81 elements. 

4. Element ID numbers must be unique with respect to all other element 
ID numbers. 

5. The field width for AREAC and THICK are 16 columns each. 
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CINGEN DATA CARD DESCRIPTION 


Card : F21 

Description: One-way fluid flow element 


Format and Example: 


F21 EID 

MID Jl J2 AREA 

F21 101 

3 12 6 11. 

Field 

Contents 

F21 

F21 signals a 2-node, one-i 

EID 

element ID 

MID 

material ID 

Jl 

upstream grid point 

02 

downstream grid point 

AREA 

cross-sectional flow area 

Remarks : 



(6A8,free field) 


1. Grid point Jl must be the upstream grid point. 

2. The conductor value is calculated as AREA x Cp x COND. 

Cp and COND are obtained from material ID. COND is the flow 
rate per unit surface area. 

3. Element ID numbers must be unique with respect to all other 
element ID numbers. 


4. The field length for AREA is 16 columns. 


CINGEN DATA CARD DESCRIPTION 


Card : C21 , C31 , or C41 
Description : Convection elements 
Format and Example : 


C21 

EID 

NFILM 

J1 

. 02 

CIRC 

(7A8,free field) 

C31 

EID 

NFILM 

01 

02 

03 


C41 

EID 

NFILM 

01 

02 

03 04 


C41 

10 

17 

1 

8 

7 2 



Field 

C21 

C3T 

C41 

EID 

NFILM 

01, 02, J3, J4 
CIRC 


Contents 

2- node convection element 

3- node convection element 

4- node convection element 
element ID 

convection property ID 

Grid points defining surface of element to which the 
convection is applied 

circumference of the element 


Remarks : 

1. Grid points must be input in a consecutive order around the element. 

2. Element ID numbers must be unique with respect to all other element 
ID numbers. 

3. An arithmetic surface node is created on the element surface to 
which the C-type element is applied except when applied to surface 
elements and 2 grid point elements. 

4. Only one surface node is created regardless to the number of C-type 
elements applied to the surface. 

5. The field length for CIRC is 16 columns. 
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CINGEN DATA CARD DESCRIPTION 


Card : R21 , R31 , or R41 

Description : Radiation elements 

Format and Example : 


R21 

EID 

NRAD 

01 

J2 

CIRC 

(7A8,free field) 

R31 

EID 

NRAD 

Jl 

J2 

J3 


R41 

EID 

NRAD 

Jl 

J2 

J3 J4 


R31 

12 

120 

7 

8 

9 



Field 

R21 

R31 

R41 

EID 

NRAO 

Jl, J2, J3, J4 
CIRC 


Contents 

2- node radiation element 

3- node radiation element 

4- node radiation element 
element ID 

radiation property ID 

grid points defining surface of element to which the 
radiation is applied 

circumference of the element 


Remarks : 

1. Grid points input in any consecutive order around the element. 

2. Element ID numbers must be unique with respect to all other element 
ID numbers. 

3. An arithmetic surface node is created on the element surface to 
which the R-type element is applied except when applied to surface 
elements and 2 grid point elements. 

4. Only one surface node is created regardless to the number of R-type 
elements applied to the surface. 

5. The field length for CIRC is 16 columns. 
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CINGEN DATA CARD DESCRIPTION 


Card: CBAR, CTRI, CQUAD, CWEDGE , CHEX 

Description : NASTRAN conduction element types for 2 , 3 , 4, 6 and 8 grid points 

Format and Example : 


CBAR 

EID 

PIDB 

01 

02 





CTRIA2 

EID 

PIDT 

01 

02 

03 




CQUAD2 

EID 

PIDQ 

01 

02 

03 

04 



CWEDGE 

EID 

MID 

01 

02 

03 

04 

05 

06 

CHEX 

EID 

MID 

01 

02 

03 

04 

05 

06 


07 

08 







CWEDGE 

7 

18 

1 

9 

7 

6 

20 

22 


Contents 

CBAR, 2-node conduction element 
CTRI, 3-node conduction element 
CQUAD, 4-node conduction element 
CWEDGE, 6 -node conduction element 
CHEX, 8 -node conduction element 
element ID 

property ID card number (PBAR card) 
property ID card number (PTRIA2 card) 
property ID card number (PQUAD2 card) 
material ID 

03, 04 grid points defining the element 
07, 08 

Remarks : 

1. Element ID numbers must be unique with respect to all other element ID 
numbers. 

2. NASTRAN and SPAR elements cannot be mixed. The type of elements 
(NASTRAN or SPAR) processed is determined by the 1 st INPUT RECORD 
(See page 3 - 13 .) 


Field 

CBAR 

CTRIA2 

CQUAD2 

CWEDGE 

CHEX 

EID 

PIDB 

PIDT 

PIDQ 

MID 

01 , 02 , 
05, 06, 
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C INGEN DATA CARD DESCRIPTION 


Card: END DATA 

De scription : This card ends the input data. 

Format and Example : 

VALUE (A8, free field) 

ENDDATA 


Field Contents 

VALUE ENDDATA to denote the last card in the input set. 

Remarks: 

1. This data card must he input. 
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The GEOMPLT program is an interactive graphics program developed, 
initially, for the display of finite element data. It has now been 
modified and expanded to include plotting a MITAS or SINDA thermal network, 
which has been created through the use of the CINGEN program (see Section 
2.2.0). In addition, an XYPLOT capability has been added for displaying 
MITAS or SINDA results in the form of Time- Temperature or Spatial - 
Temperature graphs. 

A functional flow diagram of the GEOMPLT program is illustrated in 
Figure 9. Figure 10 is a detailed flow diagram of the XYPLOT logic. It 
is recommended that the user acquaint himself with the overall logic flow 
of the program and its parameter options before entering into a graphics 
session on a terminal. However, the program will prompt the user with 
questions and options which are generally self explanatory. Best results 
are most often obtained after carefully planning the view angles, partial 
view selection and labeling options. 

Appendix B contains a user guide to the GEOMPLT program as it existed 
before the recent improvements for thermal network and XY plotting. This 
information has been included for completeness, however, the following 
sections describe the capabilities, limitations, INPUT requirements and 
OUTPUT options necessary to perform finite element, thermal network and XY 
plotting. 

4.1 CAPABILITY 

The GEOMPLT program has the following general capabilities: 

1. Graphically displays finite element models. 

2. Graphically displays thermal networks. 

3. Prepares data file for use by CQNTUR program. 

4. Edits finite element data file (CDC only). 
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If PROGRM • 7 


Plot Finite Element 
Model 


Load JPV Array with -1 


User Input: 



Write error message 
and return $XYPL0T 


Space ■ 

m 

Hist • 

a 

JPV(l)- 

Times » 
IX 
IY 
IZ 

Node ■ 
STARTT- 
ENDT ■ 
IN 


0 no spatial -temp plot 

1 spatial -temp plot 

0 no time-temp plot 

1 tlma-temp plot 

List of nodes for spatial -temp 
plot (Individually or In a string) 
spatial -temp plot time 
1 plot temp versus x-coord 

1 plot temp versus y-coord 

1 piot temp versus z-coord 

Node no. for temp-time plot 
Start time for temp-time plot 
End time for temp-time plot 
1 plot temp versus node no. 



Read unformatted file 99 
No. of nodes(MITAS FILE) 


Load JPV with nodes or string 
of nodes to be plotted in 
spatial -temp plot 


Read 99 for each node 
no, and check with 
value In NODE, 


Load arrays for X, Y, or Z coord, and node nos. 
Node nos. may be loaded Individually or In a 
string from file 29 (GINGEN written file) 


-3 













>> 



rite Error Message 


m 



>e 10, (continued) 
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''if Node No. 1 
vNODE 

« 

From file 99 ■' 
V 


Yes 

Store sequence of NODE 



Read time and tempera- 
j tures f ran file 99 


, 

< If tiro 

> END! 


No 

< ’C If time 

> STARTT 


No 

Place all temp value Into 
dummy variable. Do not 
-sav^. _ 



Figure 10* (continued) 
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5. Creates an alternate print file with diagnostic information 
for review after terminal session. 

6. Writes a data file for obtaining hard copies of geometric 
plots via Cal comp or Houston Instrument drum plotters. 

7. Graphically displays XY plots. 

Computer Configuration - Versions of the program are being maintained for 
the CDC 6600 and UNIVAC 1100 computer systems. Both versions are essentially 
the same with minor programing changes to accommodate a variable core 
concept on the UNIVAC and CDC systems. 


1 . 

Computer 

CDC 6600 

UNIVAC 1100 

2. 

Core size required 

26K min. (user re- 
compiles main to 
expand) 

22K min. (Qyn expand 
for problem size) 

3. 

Language 

FORTRAN IV 

FORTRAN IV 

4. 

System 

Scope 3.4 or NOS 

Exec 8 

5. 

Plotter 

Tektronix 4000 Series 

Same 

S. 

Card Punch 

Optional 

Optional 

7. 

Tape Drive 

None 

None 

8 . 

Mass Storage 

1-min problem size 

Same 


(Disk) 

dependent 


9. 

Run time 

221 elem. 

441 grid - 7.03 cps, 
1 plot 

14 epu sec. 

10. 

Restart 

None 

None 


Geometric Plot - The purpose of this feature is to provide a 3-D 
view of the finite element model and/or the thermal network on a 
Tektronix CRT terminal . Options may be selected through the NAMELIST 
input. Two NAMELIST's are used; Table I lists the options available under 
the $INPUT NAMELIST and Table II lists the options available under the $PVIEW 
NAMELIST. Table III lists element names which are recognized and plotted 
by the GEOMPLT program. Specific capabilities include: 


1. Graphically display finite element,CINGEN (MITAS models), and $ D AR models 

2. Display grid, element, thermal node and conductor number. 

3. Display the conductor and thermal node type with an alpha character 
preceding the ID number. 

4. Rotate any via* about three axes. 

5. Plot partial views of a finite element model containing all of one 
type element by simply naming the type in a list, plot partial 
views by providing a list of element .numbers, or a combination of 
both. 

6. Change plot scale factors for increasing or decreasing the size 
at the image plotted. 

Data Editor - The data editor provided in the CDC Version is limited 
in its capability, and is not recommended for general use. The CDC or UNIVAC 
system TEXT EDITOR is recommended. This editor can be used in a demand 
mode to correct errors found through plotting the finite element or thermal 
network data contained on files 4 and 29 respectively. 

Hard Copy - This feature is used to obtain high quality hard copies 
of the plot displayed on the screen. Appendix B» section 3.3, contains 
a detailed discussion on this option. 

CONTUR Plot Pre-Processor - This feature is used tc create two data 
files (10 and 11 which contain transformed grid coordinate data and connection 
date) for use by the CONTUR program. Appendix B, section 3.4, contains a 
detailed discussion on this option. 

Alternate Print File - This feature provides a permanent record of the 
plotting session. File 7 (CDC) is used and the user must divert the file to a 

printer if a copy is desired. Appendix B, section 35, contains a description 
of information contained on this file. 
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X-Y Plotting - This feature is used to display analytical results from 
a MITAS (or SINDA) model which was developed using the CINGEN program. An 
unformatted tabular data file must be written through execution of MITAS 
subroutine TSAVE. This subroutine copies thermal data to tape 'PLOT' which 
then must be copied to file 99 for plotting. This data file may contain 
any thermal parameter (e.g. temperature, heat flux) in tabular form as a 
function of: 

1. Time 

2. Node 

3. Spatial coordinate (X, Y, or Z) 

Options are selected thorugh the use of $XYPL0T NAMELIST input. Table IV 
lists the options available under the $XYPL0T NAMELIST. Formatted tabular 
data may be input on file 30. XY data is selected from file 30 or 99 
using the IN30 parameter in $XYPL0T NAMELIST (default is file 99). 

A typical procedure is as follows: 

1. Generate spatial coordinates (file 29) in CINGEN for use by 
plotting program. 

2. Create thermal results (file 99) with MITAS (or SINDA) for XY 
plotting. 

3. Attach file 29 and 30 or 99 for use by GEOMPLT. 

4. Begin interactive XY plotting. Note parameter PROGRAM in 
$ INPUT NAMELIST must be set equal to 7 for XY plots. 

See example on page 4-10. 
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Table I. Namelist $INPUT Parameter Definition (CDC ONLY) 


Parameter 

Description 

Default 

Other 

BAUD* 

Transmission rate from computer to 
terminal, BAUD/10 

30 

Integer 

CONDNO 

Print conductor numbers for thermal 
network 

O-NO 

1-YES 

CONTUR 

Saves XT plotter coordinates on file KC 
and connector data on file KD to be read 
by COUTUR program. , NOTE: When CONTUR = 1, 

the following default changes are used 
KTEST-" 1 
KD » 11 
MPLOT = 0 
MPLOT = 1 

0-N0 

1 or 9-YES 

END 

Print element numbers at element centroids 

O-NO 

1-YES 

GNO 

Print grid numbers at grid point 

O-NO 

1-YES 

IHUSTN 

Write hard copy plot commands* on file 25; 
request to permanently save plot on file 
26 is made aflei. plot Is complete 

O-NO 

1-YES 

RC 

File number to use when CONTUR = 1 

10 

Integer 

KD 

File number to use when LTEST or KTEST » 1 

7 

Integer 

KTEST 

«r 

Option for CONTUR element punch option 
if KTEST = 9; plotting is suppressed 

O-NO 

1 or 9-YES 

LTEST 

Option for CONTUR XT data punch options 
if LTEST <= 9i plotting is suppressed 

O-NO 

i 

1 or 9-YES 

MPLOT 

Option for multiple plots. Set MPLOT ® 0 

1-YES 

0-N0 


to terminate plotting 



NCHBDY 

NAS TRAN thermal 

elements to be plotted 

O-NO 

1-YES 

m 

HODENQ 

Print capacitance node number for thermal 
network 

0-H0 

1-YES 


MPLOT 

Number of plots 

if MPLOT = 1 

10 

Integer 

V 

KTEST 

Option to print 
on file 7 

rotated coordinate data 

S'SSS ® 3 

O-NO 

1-YES 
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Table I. (Continued) 


Parameter 

Description 

Default 

Other 

PROGRM 

Type of input data on file 4, 29, 30 or 99 

1 = Finite element (file 4) 

2 = Thermal network (CINGEN) only plotted (file 29) 

3 = CINGEN and finite element data 

plotted (file 4 and 29) 

7 = CINGEN and MITAS XY-PLOT (file 29 and 30 or 99) 

1 

1-5 

PV 

Option for partial views; if PV - 1 
Namelist $JPVLST is used to specify 
element numbers to be plotted 

0-N0 

1-YES 

SCALEX 

Option for adjusting the apparent width 
of the screen 

10 

REAL 

SCALEY 

Option for adjusting the apparent height 
of the screen 

7 

REAL 

SYMB 

Option to place a symbol at the grid point 
29 =>* 

29 

Integer 

X* 

Angle of rotation about x-axis (deg) 

34.27 

REAL 

y* 

Angle of rotation about y-axis (deg) 

23.17 

REAL 

Z* 

Angle of rotation about z-axis (deg) 

43.0 

REAL 

XAXIS 

Axis to be placed in screen x-axis 
(+ left to right) 

1 

1-6 

YAXIS 

Axis to be placed in screen y-axis 
(+ up) 

2 

1-6 

ZAXIS 

Axis to be placed in screen z-axis 
(+ out of screen) 

Structure yes may be interchanged with screen 
axes for changing view, where: 1 = x, 2 = y, 

3 - z, 4 - -x, 5 = -y, 6 = -z 

3 

1-6 

I DASH* 

Option for dashed lines in thermal network 
plot. Default is straight lines. 

0 

1 


Notes : 


1. Parameters marked by * retain input value until changed by a new input value. 

2. The PROGRAM parameter must be set equal to 7 for XY plotting. All other 
parameters in $INPUT NAMELIST are used for geometry only. See Table II for 
$XYPL0T NAMELIST parameters which may be selected for XY plotting. 
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Table II. Namelist $PVIEW Parameter Definition 


JPV(l) = List of elements to be displayed in a partial view. 

Elements to be displayed in one of the following ways: 

1} Individually 

2) A set of consecutive element numbers in increasing order 
by naming the first and last number in the set. 

3) A set of element numbers of the same type. 


EXAMPLE: JPV(l) = 1 , 5, 9, 13, -0, 18, 

6HK81 $ 

will plot element numbers 1, 5, 9, 13, 14, 15, 16, 
17, 18 and all K81 elements. 


Table III 


SUMMARY OF FINITE ELEMENTS RECOGNIZED BY GEOMPLT 


NASTRAN Type 


SPAR Type 


CBAR 

CONROD 

CROD 

CTUBE 

CTRIA1 

CTRIA2 

CTRMEM 

CTRPLT 

CQDMEM 

CQDPLT 

CQUADl 

CQUAD2 

CSHEAR 

CTWIST 

OTETRA 

PLOTEL 

CTRAPR 

CTRIAR 

CTORDR 

CWEDGE 

CHEXA1 

CHEXA2 

CIS3D8 

CIS3D20 

CHBDY 


K21 

C21 

R21 


F21 

K31 

C31 


R31 

K41 

C41 

R41 


K61 

K81 
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Table IV. $XYPLOT NAMELIST Parameters 


PARAMETERS **************** DEFINITION ****************** DEFAULT ******0THER ***** 


SPACE 

SPATIAL PLOT 

O-NO 

1-YES \ 

HIST 

TIME-HISTORY PLOT 

O-NQ 

1-YES 

IN30 

PLOTS FORMATTED DATA 

O-NO 

1-YES 

IX 

PLOT NODAL X-COORDINATES 

O-NO 

1-YES 

IY 

PLOT NODAL Y-COORDI NATES 

O-NO 

1-YES 

IZ 

PLOT NODAL Z-COORDI NATES 

O-NO 

1-YES i 

IN 

PLOT NODE NUMBERS 

O-NO 

1-YES 

NODE 

NODE NO FOR TIME-HISTORY PLOT 

BLANK 

INTEGER 

TIMES 

INSTANT OF TIME FOR SPATIAL PLOT 

0.0 

REAL 1 

STARTT 

STARTING TIME FOR TIME-HISTORY PLOT 

0.0 

REAL 

ENDT 

ENDING TIME FOR TIME-HISTORY PLOT 

0.0 

REAL 

SYMB 

SELECTS DATA POINT SYMBOLS 

1-0 

INTEGER 


***************************;Mr*************:H?******************£********************** 


JPV IS AN ARRAY CONTAINING THE NODE NOS TO BE PLOTTED. 

(-0) IS USED FOR CONSECUTIVE NUMBERS. INPUT EXAMPLE: 

JPV (1) = 1, 2, 3, 4, -0, 8, 17, 20 (PLOTS NODE NOS 1 THRU 8, 17 and 20) 


4.2 LIMITATIONS 


1. Only finite elements in Table III may be plotted. 

2. The number of elements plotted depends upon the size 
of core available. 

Core = Code + Data + Open Core 
Code + Data = 21000 d 

Open Core = 4 x NGD + 4 x NEL + JSIZE + TO + NJPV 
NGD = No. of Grid Points 
NEL = No. of Elements 

JSIZE = No. of Connections = (NEL x No. Grid PTS/Elem.) 

NJPV = (300 - NEL/ 2) if NEL < 600 
= 0 if NEL > 600 

3. Program stops if 100 plots are generated. 

4. If contour data is required for several views, a new execution 
must be initiated for each view. 

5. Up to three conductor labels can be plotted between any two 
nodes. 

6. Labeling will be overwritten if two or more nodes occupy 
the same spatial coordinates. 

7. Maximum number of time increments which may be plotted 

must not exceed LENGB/12, where LENGB = 7000. This parameter 
may be changed in the MAIN program. 


4.3 INPUT/OUTPUT DESCRIPTION 


As previously discussed, the user's response during execution of 
the GEOMPLT program is generally an answer to a question or a NAMELIST 
input. The parameters which may be selected or redefined in the 
NAMELIST arguments are described in Tables I, II, and IV. The user 
is cautioned that some parameters in the $INPUT and $XYPL0T NAMELISTS 
revert back to their original default value while others retain their 
new value until changed. All parameters have a default value which is 
listed in Table I, II, or IV, 

The Table V lists all files (input and output) which may be used 
by GEOMPLT. 
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Table V. File Allocation 



File Number 




4 

Finite element data 

Req*d 

Input 

5 

TEKTRONIX Keyboard 

Req'd 

Input 

6 

TEKTRONIX Screen 

Req'd 

Output 

7 

Alternate print file (CDC only) 

Optional 

Output 

8 

EDIT file for Finite Element Data (CDC only) Optional 

Output 

10 

Grid data for CONTUR 

Optional 

Output 

11 

Connection data for CONTUR 

Optional 

Output 

21 

Alternate print file (UNIVAC only) 

Optional 

Output 

25 

Hard copy Scratch file 

Optional 


26 

Plot Commands for Hard Copy 

Optional 

Output 

29 

CINGEN Thermal Network data (MITAS model) 

Optional 

Input 

30 

FORMATTED XY - DATA 

Optional 

Input 

99 

UNFORMATTED MITAS XY - DATA 

Optional 

Input 



GEOMPLT FILE DESCRIPTION 


FILE : 4 

Description: Finite element data cards 


Record Description 


Format 


1 

2 

3 

4 


Last 


Grid 6ID XI X2 X3 

K21 EID MID J1 J2 A 

F21 (See Section 2.1.3 "SPAR DATA 

C21 CARD DESCRIPTION", a description 

R21 of all Element output card 

K31 formats are contained under "REMARKS.") 

C31 

R31 

K41 

C41 

R41 

K61 

K81 

(K81 Continuation) 

ENDDATA 


(A4, 4X, 18, 8X, 3F8.2) 
(A3, 5X, 418, El 6. 7) 


(A7) 


REMARKS : 

1. Only SPAR Thermal Elements are shown above, however, NASTRAN element 
types shown in Table III are also permitted. 

2. Elements may be in unsorted order. 

3. If GRDSET card is used it must preceed the first GRID card. 

4. If local coordinate definition cards are used, they 
must appear before the GRID cards. 

5. All continuation cards must follow their parent. 

6. The last card image in the file must be an ENDDATA card. 

7. File 4 may contain additional records to be used by the CINGEN 
program. These will be ignored by GEOMPLT. 

8. If a CINGEN thermal network is to be displayed the thermal node 
numbers must be equal to the finite element numbers. 
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Description : Thermal Network Data Cards 


RECORD DESCRIPTION 


FORMAT 


1 TY, M, TNI, J, TN2, K, GV, X, Y, Z, A, B, 

2 TY, M, TNI, J, TN2, K, GV, X, Y, Z, A, §, 

3 TY, M, TNI, J, TN2, K, GV, X, Y, Z, A, B, 

4 TY, M, TNI, J, TN2, K, GV, X, Y, Z, A, B, 

TY, M, TNI, J, TN2, K, GV, X, Y, Z, A, B, 

TY, M, TNI, J, TN2, K, GV, X, Y, Z, A, B, 

TY, M, TNI, J, TN2, K, GV, X, Y, Z, A, B, 


LAST TY 

, M, TNl, J, TN2, K, GV, X, Y, Z 

FIELD 

CONTENTS 

TY 

ELEMENT TYPE (eg: K, F, C or R) 

M 

ELEMENT NO. 

TNI 

NODE 1 TYPE (eg: A, B, or D) 

J 

NODE 1 NO. 

TN2 

NODE 2 TYPE (eg: A, B, or D) 

K 

NODE 2 NO. 

GV 

ELEMENT CONDUCTOR VALUE 

X, Y, Z 

NODE 1 SPATIAL COORDINATES 

A, B, C, 

NODE 2 SPATIAL COORDINATES 

REMARKS: 



C 

C 

C 

C 

C 

C 

C 

c 


{3(A1 ,17) , 7F8.3) 
(3(Al,I7), 7E8.3) 
(3(A1,I7), 7F8.3) 
(3(A1 ,17) , 7F8.3) 
(3{A1,I7), 7F8.3) 
(3 (AT ,17) , 7F8.3) 
(3{A1 ,17) , 7F8.3) 
(3(A1,I7), 7F8.3) 


1. The records shown above are generated by the CINGEN program. 

2. File 29 contains thermal element types and numbet _ , thermal 


node types and numbers, element conductor values, and node coordinates. 

3. The values of J and K (thermal node numbers) must be the same as 
the finite element ID numbers (ElD'S). 


FILE; 

30 


Description: XY Data (Formatted) 


RECORD 

DESCRIPTION 

FORMAT 

1 

NM, N(1 ), N(2) , N{3) N(NM) 

(1018) 

2 

TIME (1), DATA (1), DATA (2)... DATA (NM) 

(10F8.3) 

3 

# 

TIME (2), DATA (1), DATA (2)... DATA (NM) 

(10F8.3) 

* 

LAST 

TIME (K) 9 DATA (1), DATA (2)... DATA (NM) 

(10F8.3) 


FIELD 

CONTENTS 

NM 

number of nodes 

N(i) 

i th node ID number 


REMARKS ; 

1. File 30 contains formatted data used for XY - Plotting, 

2. To select file 30 for XY data input, reset parameter IN30 = 1 
in $XYPL0T NAMELIST. 
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FILE : 99 

Description : MITAS OUTPUT (UNFORMATTED) 


RECORD 

DESCRIPTION 

FORMAT 

1 

NM, N(1 }» N(2), N(3),. N(NM) 

NONE 

2 

TIME (1), DATA (1), DATA (2),.. .DATA(NM) 

NONE 

3 

a 

TIME (2), DATA (1), DATA (2) . .DATA(NM) 

NONE 

• 

LAST 

TIME(K), DATA (1), DATA (2) . .DATA(NM) 

NONE 


FIELD 

CONTENTS 

NM 

number of nodes 

N(1) 

i th node ID number 


REMARKS: 

1. File 99 contains unformatted data used for XY - Plotting. 

2. The XY Plotting routine will read file 99 for XY data unless 
parameter IN30 is reset to IN30 = 1 in the $XYPL0T NAMELIST. 

3. See MITAS (or SINDA) Users Manual for a discussion of variables 
which may be written on file 99,and the method used. 
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Typical XY Output 

The following example assumes the GEOMPLT program has been loaded 
on the machine and all necessary data files (4 for finite element model, 
29 for CINGEN thermal network, and 99 for MITAS XY data tables) have been 
attached. After execution, the program Is conversational and the user 
must respond to requests by the program in a step - by - step manner. 

The > symbol is used to identify response by the user. 

This example displays XY plots for temperature vs. node numbers, 
temperature vs. time, temperature vs. X-coordinate, temperature vs. Y- 
coordinate, and temperature vs. Z-coordinate. Multiple plots are dis- 
played by typing a 'blank 1 (or any other character) after each plot 
appears. After 'blank' is transmitted the screen is. cleared and the 
lines from 'TYPE (C) FOR A LIST OF $XYPL0T NAMELIST OPTIONS' will appear. 
After typing an 'A* for a list of $XYPL0T NAMELIST options or a 'C' for 
no list the user must input the correct parameter values to $XYPL0T 
NAMELIST to obtain the desired type of plot. To terminate XY plotting 
the word 'STOP' must be typed and transmitted. 


THE FOLLOWING TABLE ALLOWS THE USER TO VARY THE CORE 

USAGE IN ORDER TO ALLOW PLOTTING OF LARGER MODELS 

ESTIMATED NO. OF 
ELEMENTS TO BE PLOTTED 
IN ANY PARTIAL VIEW 

0.0 

121.0 

241.0 

357.0 

471.0 

579.0 

TYPE IN ONE OF THE ABOVE FLOATING POINT NUMBERS. 

TYPE IN MAX NUMBER IF X-Y GRAPHS ARE DESIRED. 

> 579.0 

MAXIMUM NO. OF ELEMENTS THAT CAN BE PLOTTED IN A COMPLETE VIEW IS * 583 
MAXIMUM NO. OF ELEMENTS THAT CAN BE PLOTTED IN A PARTIAL VIEW IS * 579 
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TYPE IN TITLE CARD 
> (USER TYPES IN PLOT TITLE) 

TO HAVE DEFINITIONS OF 'INPUT' NAME LIST' VARIABLES DISPLAYED, TYPE IN C 
OTHERWISE, TYPE IN ANY OTHER CHARACTER 


> C 


GNO = 1 TO PRINT GRID NUMBERS. DEFAULT = 0 
ENO = 1 TO PRINT ELEMENT NUMBERS. DEFAULT = 0 
MPLOT = 1 FOR MULTIPLE PLOTS. DEFAULT = 1 
NPLOT = NUMBER OF PLOTS IF MPLOT * 1. DEFAULT = 10 
PV = 1 FOR PARTIAL VIEW OPTION. DEFAULT = 0 
CONTUR = 1 SAVES XY PLOTTER COORDINATES ON FILE KC AND 
CONNECTOR DATA ON FILE KD TO BE READ BY CONTOUR 
PLOTTING PROGRAM. DEFAULT = 0 
NOTE-UP TO 2000 GRIDS ARE ALLOWED PER CONTOUR PLOT. 

THE FOLLOWING DEFAULT CHANGES ARE USED WHEN CONTUR = 1 
KTEST - 1 KD = 1 1 

NOTE-WHEN CONTUR = 1 , MULTIPLE PLOTS ARE NOT ALLOWED 
MPLOT IS SET = 0, NPLOT IS SET = 1 
CONTUR = 9 FOR CONTUR PLOTS ONLY - SAME AS CONTUR = 1 
EXCEPT GEOMETRY PLOTS ARE SUPPRESSED. 

X=ANGLD = ANGLE OF ROTATION ABOUT BASIC X- AXIS 
DEFAULT = 34.27 

Y = ELEV = ANGLE OF ROTATION ABOUT BASIC Y - AXIS 
DEFAULT =23.17 

Z =ROTX = ANGLE OF ROTATION ABOUT BASIC Z-AXIS 
DEFAULT =0.0 

IDENTIFY TYPE OF DATA TO BE PLOTTED, OPTIONS ARE: 

PROGRM = 1 FOR FINITE ELEMENT DATA 
PROGRM = 2 FOR THERMAL NETWORK (CINGEN) DATA 
PROGRM = 3 FOR FINITE ELEMENT + CINGEN DATA 

PROGRM * 7 FOR CINGEN t MITAS XY PLOT 

BAUD = BAUD RATE (DEFAULT = 30) 

XAXIS = AXIS TO BE PREPOSITIONED IN BASIC X DIRECTION , (+ UP THE PAGE) 

YAXIS « AXIS TO BE PREPOSITIONED IN BASIC Y DIRECTION , (+ RIGHT TO LEFT) 

ZAXIS = AXIS TO BE PREPOSITIONED IN BASIC Z DIRECTION , (+ OUT OF PAGE) 

NOTE - LOOKING TOWARD ORIGIN FROM POSITIVE END OF EACH AXIS 

CLOCKWISE ROTATIONS ARE POSITIVE 
NOTE - CODES FOR PREPOSITIONING ARE 
XX = 1, YY = 2, ZZ=3, MX =4, MY = 5, MZ = 6 

DEFAULT AXIS VALUES ARE XAXIS = 1, YAXIS = 2, ZAXIS = 3 

**TYPE IN A CHARACTER TO CONTINUE** 
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> c 

TYPE IN NAMELIST/$INPUT/ 

NQTE-THE FIRST COLUMN IS IGNORED IN NAMELIST INPUT 

THE LAST CHARACTER MUST BE A $, NOT ON A SEPARATE LINE, BUT AFTER 
LAST VARIABLE. 

> $INPUT PROGRAM = 7 $ (EXAMPLE) 

TYPE (C) FOR A LIST OF ‘XYPLQT 1 NAMELIST OPTIONS 
TYPE ANY OTHER CHARACTER FOR NO LIST 


> C 


PARAMETERS 

*********** DEFINITION ************** 

• DEFAULT 

k* OTHER ******* 

SPACE 

SPATIAL PLOT 

O-NO 

1-YES 

HIST 

TIME-HISTORY PLOT 

O-Nu 

1-YES 

IN30 

PLOTS FORMATTED DATA 

O-NO 

1-YES 

IX 

PLOT NODAL X-COORDINATES 

O-NO 

1-YES 

IV 

PLOT NODAL Y-COORDINATES 

O-NO 

1-YES 

IZ 

PLOT NODAL Z-COORDI NATES 

O-NO 

1-YUS 

IN 

PLOT NODE NUMBERS 

O-NO 

1-YES 

NODE 

NODE NO FOR TIME-HISTORY PLOT 

BLANK 

INTEGER 

TIMES 

INSTANT OF TIME FOR SPATIAL PLOT 

0.0 

REAL 

STARTT 

STARTING TIME FOR TIME-HISTORY PLOT 

0.0 

REAL 

ENDT 

ENDING TIME FOR TIME-HISTORY PLOT 

0.0 

REAL 

SYMB 

SELECTS DATA POINT SYMBOLS 

1-0 

INTEGER 


****************************************************;* **\. ' ************** 


JPV IS AN ARRAY CONTAINING THE NODE NOS TO BE PLOTTED. 

(-0) IS USED FOR CONSECUTIVE NUMBERS, INPUT EXAMPLE: 

JPC (I) = 1, 2, 3, 4, -0, 8, 17, 20 (PLOTS NODE NOS 1 THRU 8, 17 and 20) 

**TYPE ($XYPLOT (OPTIONS) $ ) 

** NOTE - 1ST CHARACTER IS IGNORED FOR NAMELIST INPUT 
** (,) AS LAST CHARACTER TO CONTINUE ON NEXT LINE 

** ($) AS LAST CHARACTER TERMINATES NAMELIST INPUT. 

> $XYPLOT SPACE = 1, TIMES = 4.0, IN = 1 , OPV (1) = 5001, -0, 5004 $ (EXAMPLE) 
AFTER PLOT IS COMPLETE TYPE IN A BLANK TO CONTINUE PLOTTING 

TYPE CTO CONTINUE 

> C 

(PLOT APPEARS) SEE FIG. 11 
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5001.0 5001.5 5003.0 5003.5 5003.0 5003.5 

NODE NUMBER 


TEST XVPLOT 


5004.8 


Figure 11, Typical XY Plot Generated fey GEOMPLT 




TYPE (C) FOR A LIST OF $XYPLOT NAMELIST OPTIONS 
TYPE ANY OTHER CHARACTER FOR NO LIST 

> A 

**TYPE ($XYPLOT (OPTIONS) $ ) 

** NOTE - 1ST CHARACTER IS IGNORED FOR NAMELIST INPUT 
** (,) AS LAST CHARACTER TO CONTINUE ON NEXT LINE 

** ($) AS LAST CHARACTER TERMINATES NAMELIST INPUT. 

> $XYPLOT HIST = 1, NODE = 5002, STARTT = 0.0, ENDT = 10. 

AFTER PLOT IS COMPLETE TYPE IN A BLANK TO CONTINUE 
TYPE C TO CONTINUE 

> C 

(PLOT APPEARS) SEE FIG. Ha. 


(EXAMPLE) 
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TYPE (C) FOR A LIST OF $XYPLOT NAMELIST OPTIONS 
TYPE ANY OTHER CHARACTER FOR NO LIST 

> A 

**TYPE ($XYPLOT (OPTIONS) $ ) 

** NOTE - 1ST CHARACTER IS IGNORED FOR NAMELIST INPUT 
** (,) AS LAST CHARACTER TO CONTINUE ON NEXT LINE 

** ($) AS LAST CHARACTER TERMINATES NAMELIST INPUT. 

> $XYPLOT SPACE = 1, TIMES ** 4.0, IX « 1, JPV (1 ) • 5001, -0 5004 $ (EXAMPLE) 
AFTER PLOT IS COMPLETE TYPE IN A BLANK TO CONTINUE 

TYPE C TO CONTINUE 

> C 


(PLOT APPEARS) SEE FIG. 11b. 


mm x womimiE 


TEUP US X-Ce@R9IftATE XVPLOT EXAMPLE 

figure lib. Typical XY Plot Generated by GE0MPL1 


TYPE (C) FOR A LIST OF $XYPLOT NAMELIST OPTIONS 
TYPE ANY OTHER CHARACTER FOR NO LIST 

> A 

**TYPE ($XYPLOT (OPTIONS) $ ) 

** NOTE - 1ST CHARACTER IS IGNORED FOR NAMELIST INPUT 
** (,) AS LAST CHARACTER TO CONTINUE ON NEXT LINE 

** ($) AS LAST CHARACTER TERMINATES NAMELIST INPUT. 

> $XYPLOT SPACE = 1, TIMES = 4.0, IY = 1, JPV (1) = 5001 , -0 5004 $ (EXAMPLE) 
AFTER PLOT IS COMPLETE TYPE IN A BLANK TO CONTINUE 

TYPE C TO CONTINUE 

> C 

(PLOT APPEARS) SEE FIG. 11c. 
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TYPE (C) FOR A LIST OF $XYPLOT NAMELIST OPTIONS 
TYPE ANY OTHER CHARACTER FOR NO LIST 

> 

**TYPE ($XYPLOT (OPTIONS) $ ) 

** NOTE - 1ST CHARACTER IS IGNORED FOR NAMELIST INPUT 
** (,) AS LAST CHARACTER TO CONTINUE ON NEXT LINE 

** ($) AS LAST CHARACTER TERMINATES NAMELIST INPUT. 

> $XYPLOT SPACE = 1, TIMES = 4.0, IZ = 1, JPV (1) - 5001,-0 5004 $ (EXAMPLE) 
AFTER PLOT IS COMPLETE TYPE IN A BLANK TO CONTINUE 

TYPE C TO CONTINUE 

> C 

(PLOT APPEARS) SEE FIG. lid. 


iT'V“-ij xf "¥‘ r V~ 1 u 'Jl "J Tu'Tf j 8 u u U jj u U u u y u 

) ^Qo@ ©oil) o © - &o'0 


zmm. e mmmKm. 


trn^n sot 


wlwt mmmic, m& 


Figure lid. Typical XY Plot Generated by GEOMPLT 
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Section 5 


SAMPLE PROBLEM 

This section presents a thermal problem which was modeled using 
SPAR finite elements. The problem (Figure 12) consists of conduction, 
convection, radiation, and fluid flow. This problem was created only to 
illustrate program capability and has no real physical identity. Data 
input and output for the three programs used in automated data generation 
are presented. 

5.1 SPAR MODELING 

The physical system was represented by SPAR finite elements. Every 
type of element available to the user is included in the problem. Figure 13 
is a listing of the SPAR input describing the model. Figure 14 shows SPAR 
output to file 4. Grid points and element card are contained on this file. 

5.2 CINGEN MODELING 

Before data file 4 can be used by CINGEN to produce a MITAS data deck, 
property and end cards must be added. Figure 15 shows file 4 after these 
cards have been added. These cards are preceeded by the word "New" in 
the listing. Figure 16 is the MITAS deck output by CINGEN. Figure 17 
is a listing of file 29 which is used by GEOMPLT. 

5.3 MODEL CHECK-OUT 

Check-out of the thermal network created by SPAR/CINGEN can be accomplished 
by the use of GEOMPLT. Several types of plots can be generated to assure 
correct modeling. Some typical plots are presented. The first page of each 
plot gives the required input to generate the plot that follows it. 
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Figure 13 


SPAR Input For Sample Problem 
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Figure 13 (tontinued) 
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Figure 14. SPAR Output to File 4 for Sample Problem 
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Figure 14. (continued) 


000113 

ODD 

RJ1 

46 

4 

41 

000114 

000 

K21 

47 

2 

46 

000115 

000 

K22 

46 

2 

47 

000116 

000 

K21 

49 

2 

4B 

000117 

000 

K2X 

50 

5 

50 

ODDI IB 

000 

K21 

51 

5 

51 

000119 

000 

K21 

52 

5 

52 

000120 

000 

K21 

53 

5 

53 

000121 

000 

IC21 

54 

5 

54 

000122 

000 

K21 

55 

5 

55 

000123 

000 

F21 

56 

3 

57 

000124 

000 

F21 

57 

3 

50 

000125 

000 

C21 

58 

5 

47 

000126 

ODD 

C21 

59 

6 

43 

000127 

GOO 

R21 

60 

5 

47 

000126 

000 

R21 

61 

6 

46 


tn 

'•J 


43 45 

47 5.6900000-01 

48 5.8900000-01 

49 5.8900000-01 

51 0.0000000 

52 0.0000000 

53 0.0000000 

54 0.0000000 

55 0.0000000 

56 0.0000000 

56 1.9600000-01 

59 1.9600000-01 

46 1.5700000+00 

49 1.5700000+00 

48 3.1400000+00 

49 3.1400000+00 



Figure 14. (continued) 
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Figure 17. Data File 29 for Sample Problem as Created by CINGEN 
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Figure 17. (continued) 
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** TYPE TITLE CARD <80 CHARACTERS) 

> SAMPLE FINITE ELEMENT DATA 
** TYPE <C) FOR A LIST OF NAMELIST OPTIONS 
** OR TYPE ANY OTHER CHARACTER FOR NO LIST 
>A 

*# TYPE ( S INPUT (OPTIONS) S) 

NOTE- 1ST COLUMN IS IGNORED FOR NAMELIST INPUT* 

** < # ) AS LAST CHARACTER T© CONTINUE ON NEXT LINE* 

mi m LAST CHARACTER TERMINATES NAMELIST INPUT 
> S1NPUT X»-S0«*V«15«*Z'»e.*SCALEX-ia« # SCALEY<'"lS« 9 


5-17 



— 


•I V 

r 


* 

♦ * TVPE TITLE CARD (80 CHARACTERS) 

CINGEN + FINITE ELEMENT DATA, K81 K61 
** TVPE ( •INPUT (OPTIONS) *) 

** NOTE- 1ST COLUMN IS IGNORED FOR NAMELIST INPUT, 

** ( , ) AS LAST CHARACTER TO CONTINUE ON NEXT LINE, 

** (•) AS LAST CHARACTER TERMINATES NAMELIST INPUT 

•INPUT NODENO - 1 , CONDNO - 1 , PU - 1 • 

** TVPE ( •PUIEU JPU(l)- , (ELEMENT NUMBERS) •) 

** TO DEFINE PARTIAL UIEU 



** NOTE - 1ST COLUMN IS IGNORED FOR NAMELIST INPUT, 

** ( , ) AS LAST CHARACTER TO CONTINUE ON NEXT LINE, 

** (•) AS LAST CHARACTER TERMINATES NAMELIST INPUT 

** (-0) IS USED FOR CONSECUTIUE NUMBERS, 

** (6H( 6 CHAR ELEMENT NAME)) TO PLOT ALL OF ONE 

** TVPE ELEMENT 

** EXAMPLE JPU( 1 )“1,2,3, -0,8,6HCBAR • 

** PLOTS ELEMENTS 1,3, 3,4,5, 6,7,8 AND ALL CBARS 

•PUIEU JPU(1 )-6HK81 ,6HK61 • 
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X 

*t TYPE TITLE CARD (80 CHARACTERS) 

C INGEN + FINITE ELEMENT DATA K41 + K31 
XX TYPE ( SINPUT (OPTIONS) S) 

XX NOTE- 1ST COLUMN IS IGNORED FOR NAMELIST INPUT, 

** ( , ) AS LAST CHARACTER TO CONTINUE ON NEXT LINE, 

** <S) AS LAST CHARACTER TERMINATES NAMELIST INPUT 

SINPUT CONDNO • i , NODENO ■ 1 , PU • 1 9 
XX TYPE ( SPUIEU JPU( 1 )• , (ELEMENT NUMBERS) S) 

XX TO DEFINE PARTIAL UIEU 

XX NOTE - 1ST COLUMN IS IGNORED FOR NAMELIST INPUT, 

XX < , ) AS LAST CHARACTER TO CONTINUE ON NEXT LINE, 

X* (S) AS LAST CHARACTER TERMINATES NAMELIST INPUT 

XX (-0) IS USED FOR CONSECUTIUE NUMBERS, 

XX (6HC6 CH«* ELEMENT NAME)) TO PLOT ALL OF ONE 

XX TYPE ELEMENT 

XX EXAMPLE JPU(1)-1,2,3,-0,8,6HCBAR 9 

XX PLOTS ELEMENTS 1,2,3,4,5,6,7,8 AND ALL CBARS 

SPUIEU JPU(1 )-6HK41 ,6HK31 9 
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C INGEM + FINITE ELEMENT DATA K41 + K31 


* 

XX TVPE TITLE CARD (80 CHARACTERS) 

>C INGEN THERMAL DATA 

XX TVPE CC) FOR A LIST OF NAMELIST OPTIONS 
XX OR TVPE ANV OTHER CHARACTER FOR NO LIST 
>A 

XX TVPE ( • INPUT (OPTIONS) S) 

XX NOTE- 1ST COLUMN IS IGNORED FOR NAMELIST INPUT, 

XX (,) AS LAST CHARACTER TO CONTINUE ON NEXT LINE, 

XX <*) AS LAST CHARACTER TERMINATES NAMELIST INPUT 

> • 

0 

•INPUT X--20o ,V-lS.,Z-e.,SCALEX"13*,SCALEV-15. • 
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XX TVPE TITLE CARD (80 CHARACTERS) 

>CINGEN THERMAL DATA UITH NODE NO. 

XX TVPE (C) FOR A LIST OF NAMELIST OPTIONS 
** OR TVPE ANV OTHER CHARACTER FOR NO LIST 

>A 

XX TVPE < SINPUT (OPTIONS) •) 

X* NOTE- 1ST COLUMN IS IGNORED FOR NAMELIST INPUT, 

XX (,) AS LAS*" CHARACTER TO CONTINUE ON NEXT LINE, 

XX (S> AS LAST CHARACTER TERMINATES NAMELIST INPUT 

> SINPUT SCALEX-13.,SCALEY-15.,X— 20. , V-iS. ,Z-0. ,NODENO-l S 
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Ci'NGEM THERNAL DATA UITH NODE NO* 


matt 




ORIGINAL} PAGE IS 
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cn 
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no 

cr* 


♦t TYPE TITLE CARD <80 CHARACTERS) 

✓ C INGEN THERMAL DATA UITH CONDUCTOR NOS. 

*X TYPE <C) FOR A LIST OF NAMELIST OPTIONS 
** OR TYPE ANY OTHER CHARACTER FOR NO LIST 

>A 


** TYPE < SINPUT (OPTIONS) S> 

X* NOTE- 1ST COLUMN IS IGNORED FOR NAMELIST INPUT, 

XX ( , ) AS LAST CHARACTER TO CONTINUE ON NEXT LINE, 

XX <S) AS LAST CHARACTER TERMINATES NAMELIST INPUT 

> SINPUT CgNDNO-l,X— 20. , Y-1S. ,Z-0. ,SCALE»13. , SCALEY-1S. S 


o 



t 


$ 


o 


0 



C INGEN THERMAL DATA WITH CONDUCTOR NOS. 


* 

** TYPE TITLE CARD (80 CHARACTERS) 

FINITE ELEMENTS C41,C31,C21 UITH NODE N0. 

** TVPE ( •INPUT (OPTIONS) • ) 

X* NOTE- 1ST COLUMN IS IGNORED FOR NAMELIST INPUT, 

** ( , ) AS LAST CHARACTER TO CONTINUE ON NEXT LINE, 

XX * (•> AS LAST CHARACTER TERMINATES NAMELIST INPUT 

•INPUT PU-i,GNO-l • 

XX TYPE C •PUIEU JPUC1 )■ , (ELEMENT NUMBERS) •) 

** TO DEFINE PARTIAL VIEW 

** NOTE - 1ST COLUMN IS IGNORED FOR NAMELIST INPUT, 

** <*) AS LAST CHARACTER TO CONTINUE ON NEXT LINE, 

XX (•> AS LAST CHARACTER TERMINATES NAMELIST INPUT 

** (-0) IS USED FOR CONSECUTIVE NUMBERS, 

** <6H<6 CHAR ELEMENT NAMED TO PLOT ALL OF ONE 

XX TVPE ELEMENT 

** EXAMPLE JPU(1)«1,2,3,-8,8,6HCBAR • 

XX PLOTS ELEMENTS 1,2,3,4,5,6,7,8 AND ALL CBARS 

•PUIEU JPU(1)-6HC41 ,6HC31 ,6HC21 • 

© 



5-30 


V 




* 

tt TYPE TITLE CARD (83 CHARACTERS) 

FINITE ELEMENTS R41,R31,R21 WITH NODE NO. 

** TYPE ( • INPUT (OPTIONS) •> 

** NOTE- 1ST COLUMN IS IGNORED FOR NAMELIST INPUT, 

** (,) AS LAST CHARACTER TO CONTINUE ON NEXT LINE, 

** <«) AS LAST CHARACTER TERMINATES NAMELIST INPUT 

•INPUT GNO-l,PU-l • 

** TYPE ( SPUIEU JPVC1)- , (ELEMENT NUMBERS) •) 

** TO DEFINE PARTIAL UIEU 

** NOTE - 1ST COLUMN IS IGNORED FOR NAMELIST INPUT, 
tt < , ) AS LAST CHARACTER TO CONTINUE ON NEXT LINE, 

tt (•) AS LAST CHARACTER TERMINATES NAMELIST INPUT 

** (-0) IS USED FOR CONSECUTIUE NUMBERS, 

tt <6H<6 CHAR ELEMENT NAPE)) TO PLOT ALL OF ONE 

** TYPE ELEMENT 

** EXAMPLE JPU(1 )- 1,2,3, -0,8, 6HCBAR • 

tt PLOTS ELEMENTS 1,2,3,4,5,6,?,8 AND ALL CBARS 

•PUIEU JPU(1)-6HR41 , 6HR31 ,6HR21 • 
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Section 3 


STRUCTURE DEFINITION 


Four programs - TAB, ELD. E, and EKS — are used to generate and store in 
the data complex data sets that define the finite element model of the structure. 

TAB and ELD are used to generate the basic definition of the structure. 
Subprocessors within TAB translate card input data into tables of joint locations, 
material constants, section properties, etc. Subprocessors within ELD translate 
card input data into data tables that define individual finite elements of various 
types . 


Using the data produced in TAB and ELD, programs E and EKS generate 
an array of data sets, collectively known as the E-state, that contain a complete 
description of every element in the structure, including details of element geometry, 
intrinsic stiffness matrices, etc. 

Except where specifically noted to the contrary, all data sets produced by 
TAB, ELD, E, and EKS should be retained in library X. 


3-1 


£ 

A- 2 


3. 1 TAB - BASIC TABLE INPUT 

Function. TAB contains an array of sub-processors which 
generate tables of material constants, section properties, joint 
locations, etc., and various other data sets comprising a substantial 
portion of the definition of the structure. Each of these sub-processors 
is identified either by a multi-word (long form) name, such as JOINT 
LOCATIONS, or by a short name, e. g. , JLOC. Each sub-processor 
generates a data set having the same name as the sub-processor. 
Sub-processor names and functions are summarized in Table TAB- 1. 

TAB may be used to either (1) create new data sets, or (2) 
update existing data sets by replacing individual entries in them. The 
update mode is commonly used in Demand (teletype) interactive operation. 
Only one direct access library, SPAR logical 1, is usable in conjunction 
with TAB. When beginning a new problem , the first data card following 
@XQT TAB (and any RESET commands) must be 

START j, m j , m^, 

In the above, j = the total number of joints in the structure. It is not 
harmful to have some unused joints (i. e. joints connected to no elements), 
for convenience in interpreting the output. This should not be carried to 
extremes, however, since it wastes core storage. It Bhould be noted that 
it is not necessary to set aside large blocks of unused joint numbers in 
areas where you expect to later augment the model. The JSEQ sub- 
processor and an array of data modifier statements allow models to be 
extended without extensive repunching of existing data cards. 


TAii 


The parameters m^, m^, - - - on the START card list joint 
motion components which are identically zero for all joints. Those 
components are relative to the joint reference frames. For m=l,2, or 3 a 
direction m displacement is indicated. For m=4, 5, or 6, a direction 
m-3 rotation is indicated. 

For example, a 1000 joint space truss would begin with 
START 1000, 4,5,6 

The above indicates that all three rotation components are zero at all 
joints. As another example, a 500 joint plane frame would begin with 

START 500 3, 4, 5 

The above assumes that the the direction- 3 displacements are normal 
to the plane of the frame, and the frame deforms in-plane only. 

SPAR uses the m’s given on the START card to determine 
submatrix storage block size (see Ref.l, Sec. 4). Therefore the START 
card (rather than the CON sub-processor) should be used to suppress 
joint motions in cases such as those illustrated above. 

Following the START card, or directly following @XQT TAB if 
Library 1 already exists, the TAB card input stream is as illustrated 
below. 

Procj 
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data cards rend by ProCj 


p roc 2 


data cards read by Proc^ 


etc. 

In the above, Proc. represents a data card containing either the short 
or long-form name of a sub-processor to be executed? 6ub-prpces*or$ 
may be execoted in any order, subject to the following restriction; if 
(he data read by one sub-processor refers to another TAEUgener^fed 
table, the other table must already ejfist- It is always safe to execute 
sub-processors in the order they appear in Table TAB-1. 

TAB input card rules 

In addition to the rules generally applying to SPAR free -field 
input, the following rules apply specifically to TAB. 

(1) Trailing items omitted on data cards are assumed to b® 
zero, except when indicated otherwise in discnssione of 
specific sub-processors. 

(2) Input should be in real (fixed on floating point) foripat, exoept 
for integer items such as entry numbers, jointnumbers, table entry 
pointers, and control variables. 


TAB 


LU 


(3) Within the data card stream being read by any sub-processor, 
the following commands may be injected, 

FORMAT = j 
MOD = m 

NREF = n 

In the above , 

O j identifies which of several admis Bible 
formats is to apply to subsequent cards. 

Details are given in discussions of individual 
sub- processors. 

O m is added to the table entry indicator 

(usually called k), in all processors other than TEXT, 

JLOC, RMASS, JREF,, CON and JSEQ. In J3LOC, RMASS, 
JREF, and CON, m is added to joint numbero 
indicated on subsequent input cards, 

O In processors JLOC, JREF, and MREF, reference 
frame n applies to data on subsequent cards. 

These three commands may be used repeatedly in the card 
stream of any sub-processor. Upon beginning execution of a 
new sub-processor, these parameters are internally reset to 
their default values, FORMAT=l, MOD=0, and NREF=1. 

(A) Each input record, including continuation cards, if any, 
must not contain more than thirty words . 

3. 1-4 


A-6 


TAB 


Correcting existing data sets 

To enter the "update" mode of operation, the following command is 
injected in the input stream. 

UPDATED 

To leave the update mode, the command is 
UPDATED 

UPDATE commands should immediately precede sub-processor execution 
commands. When operating in the update mode, the output data set produced 
in the current execution is identical to that produced in the preceding execution, 
except for entries defined by the card input of the current execution. The only 
sub-processors which cannot be operated in the update mode are TEXT and 
JSEQ. 


As an example, suppose the location of joint 1742 is found to be in 
error. The JLOC data set could be repaired by the following card sequence. 

UPDATE=1 

JLOC 

1742, 947.62, 1841.9 23.487$ 
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Library Title 

The following statement will cause an identifying title to be embedded in 
library 1 : 

TITLE* GO-character alphanumeric title 

The title, which resides in a data set named NDAL , is displayed At the 

beginning of printouts of library Tables of Contents (e.g. , via DCU/TOt^),, 

RESET Controls 

No special RESET controls are available in TAB. 

Core Requirements 

Working core requirements will not significantly exceed the larger of 
(1) 13 times the number of joints, or (2) the longest table generated (see ; - } 

Table TAB-1). 

Code Release Data 

Level 9 (UNI VAC, CDC) July 1975, coded by W. D. Whetstone. 

i 

' <£> 
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Sho-rfc 

Marne 

T^bto ... 
Length ' 

■Long fo'MTl 'NaijRflMf, 

‘Function 

TEXT 

Is <tt 

TEXT. Creates a .data 'set ^obtaining 
calpha nu m eric text documenting the analysis* 

MATO 

•0 

C 0 ii 

MATERIAL CONST ANTS., Entries define 
linear material constants E, t»,, etc,. 

NSW 


DISTRIBUTED WEIGHT., Entries define 
iionstructural distributed weight parameter s 
.(weight/ length or weight/ area).. 

A&TREF 

IE 

ALTERNATE* REFERENCE. FRAMES. Entries 
describe reference frames selected by the 
Analyst for •convenience in defining joint locations, 

etcr. 

JL©C 

3 j 

joint Locations. 

&&£f 

3 

JOINT REFERENCE FRAMES. Entries define 
the jofientation o£ reference frames associated 
with joints,, used m defining Constraint, applied 
loading.,, etc,.. 

MR^F 

5 QTi 

BEAM 'ORIENTATION. Entries define beam 

c.r,(s«s^ ’Se etion orientation . 

BRL 

8 31 

BEAM RIGID LINKS, Entries define rigid 
links offsetting the end points of elastic 
jLno;de elements from, the joints they connect* 

BA 

31 n 

KU .SECTION PROPERTIES. 

BB 

31 zx 

BEAM S 6 , x r 6 i Entries define elastic 
foharaeteristj.es ©.£ type -E 22 and E25 elements-, 

BC 

d> n 

E23 SECTION .PROPERTIES* 

BD 

8 n 

£24. SECTION’ PROPERTIES. 

SA 3ae gecB. 

SHELL SECTION PROPERTIES. 

SB 

4 m 

PANEL SECTION PROPERTIES. 

CON 

h 6j 

CONSTRAINT DEFINITION. 

JSEQ 

j 

JOINT ELIMINATION SEQUENCE, 

RMASS 

% 

RIGJD MASSES. 


*n is the number ©f bntrioo ill the ‘table.,. and 3> ts the (Jntal3iumbei°:of joints in the 
structure. 


3- 1-t 


Ml TAB/ 

TEXT 

MATC 

NSW 

5.1.1 TEXT 

The TEXT subprocessor gives the analyst a means of embedding in the output 
library a date set containing alphanumeric text descriptive of the analysis 
being performed. Each card has a 4/8 punch in column 1, followed by a 60- 
character alphanumeric string. The contents of TEXT data sets may be 
printed using Che DCU/PRINT command. 


a 


r 


a„ 


3.1.2 MATERIAL CONSTANTS (MATC) 

MATC generates a tablo of material constants, 
defining the k-Ot entry in the table is k, E, v, p, 

E *= Modulus of elasticity 
V <= Poisson's Ratio 
p ** Weight per unit volume 

Thermal expansion coefficient, direction x 
thermal expansion ceefficient, direction y 


The data sequence on the card 


6, where 


1 


OLr 

i 

0 


Angle between axes of element reference frame (x,y) and (x,y). 
Element reference frame orientation is discussed in Section 3.2. 



E- = a, 
x 1 


£- » a 

y 2 

Y = 0 


temperature 

temperature 

originae is 

of poor ot1nUTY 


If G ia omitted, the program sets 0 = 0. If a_ is omitted, the program 
aeto otp => (isotropic material). 0 must be given in degrees . 


Material constants associated with specific elements are established by 
inputs to processor ELD, by pointing to entries in the MATC table. 


3.1.4 DISTRIBUTED WEIGHT (NSW) 

A table of rum-structural distributed weight parameters is defined. The data 
sequence for the input card defining the k-th entry in the table is k, W, where 

for Z'-node elements, W = weight per unit length, and 
for 3 or 4-node elements, W => weight per unit area. 

Non-stfuctural weight attached to specific elements is defined in processor 
ELD by pointing to entries in the NSW table. 
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3. 1.4 ALTERNATE REFERENCE FRAMES (ALTREF) 

In addition to the global reference frame, the analyst may find 
it convenient to define additional reference frames. These frames have 
several uses, including the following: 

(1) Joint locations may be defined in any frame the analyst 
finds most convenient (see JLOC). 

(2) Joint reference frame orientations may be defined via 
the alternate frames (see JREF). 

Each frame is uniquely identified by a positive integer. The 
global frame is always frame 1 ; accordingly the analyst is free to define 
only frames 2, 3, - - The order of data on the input card defining 
frame k is 

k, ijfSj, *2’ a 2* V a 3' X l' X 2* X 3. 

In the above Xj, x^, and x^ are position coordinates, relative to the 
global frame, of the origin of frame k. The x's need not be given if 
only the orientation of frame k is of significance, which often is the case. 

The parameters i ^ , a^, a^ indicate ordered rotations 

defining axis orientations. Two formats are provided. In both eases, we 
begin with a local frame parallel to the global frame. After the ordered 
rotations are completed, the local frame is parallel to frame k. Each of 
the i's may be 1, 2, qr 3, in any order. The a's are angles in degrees . 

If FORMAT = 1 (default) , the sequence is: (1) rotate the local 
frame aj degrees about local axis i^, then (2) from the new position, 

3. 1.4-1 
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rotate the local frame a ^ degrees about axis i^,, then (3) from the resulting 
position, rotate the local frame degrees about axis If FORMAT=2, 
the i'o and a's indicate rotation of the global frame relative to frame k. 

In the following example of an input card defining frame 27, it 
is assumed that FORMAT = 1 . 


27 3,30. 1,10. $ 



First rotation, 

G, = Global axes 
i 



K. = axes of Frame 27. 


TAB/ 

JLOC 


3.1.5 JOINT LOCATIONS (JLOC) 

JLOC produces a table containing the position coordinates of 

the joints. The data sequence on input cards is as follows: 

k, xa^, xa^, xb^, xb^, xb^, ni, ijump, nj 

If nj is given, a second card must appear, 

jjump, xcj, xc,, xc 3 » xdj, xd£, xd^ 

There are three possible interpretations of the above: 

(1) If only k, xaj, xa^ and xa^ are given, the xa's are coordinates 
of joint k. 

(2) If k, xa — - - xb^, ni, and ijump are given, the xa's and 
xb's are coordinates of points A and B terminating a string of 
ni equally-spaced joints, 
k, k+ijump, k+2(ijump) - - - - 


Point B 


Point A 
k 


k+ 

2 (ijump) 


i 

* 

k+ 

(ni-i) (ijump) 


k+ 

ijump 

The default value of ijump is 1 . 


(3) If nj is given, a linearly interpolated two-dimensional mesh of 
(ni)(nj) joints is defined, as indicated below. 

-V D 





C 

-f 

+ 


f 


k+ 

2(jjump) 

~~ + 

+ 

+ 



k+ 


Jr-f 


■f ijump+ 

+ 

jjump 


jjump 


A 

+ 

k 

+ 

k+ 

ijump 

+ 

k+ 

2 (ijump) 


4 
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TAB/ 

JLOC 


J"V. 



t 




Although the output table generated by JLOC is in rectangular 
coordinates relative to the glpbal frame, coordinate data appearing on 
the input cards may be either rectangular or cylindrical, and may be 
relative to any frame already defined in A.LTREF. The associated 
control cards are summarized below. 

Control Meaning 

NREF= n Coordinate data on subsequent cards is relative 

to frame n (until another NREF command is 
encountered). 

FORMAT = 2 Subsequent data is in cylindrical coordinates, 
relative either to frame 1 (global) or to any 
other frame selected by an NREF command. 

The convention is shown below. 

1 

Jtj = radial distance from 3 -axis . 

x_ s angle, in degrees , from 
the 1-3 plane . 

x 2 = linear distance from the 
1-2 plane. e 

FORMAT =1 Switch back to rectangular coordinates. 

Switching among frames, and between rectangular and cylindrical 
coordinates is unrestricted. 

If cylindrical coordinates are used in connection with mesh generation, 
interpolation is performed before transformation to rectangular coordinates; 
so that regular meshes on circles, cylinders and cones are readily generated. 
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A-14 



LI 1 


ELD 


3. 2 ELD - ELEMENT DEFINITION PROCESSOR 

Function - ELD translates element definition data from input cards 
into data sets which are readily usable by other SPAR processors. Elements 
may be defined singly or . through a variety of network generators, or com- 
binations thereof. An element is defined by specifying (1) the joints to which 
it is connected, and (2) pointers to applicable entries in tables of section 
properties, material constants, etc. The tables of section properties, etc., 
are usually generated by the TAB processor. As the ELD input is processed, 
checks are performed to detect errors such as references to nonexistent 
table entries, joint numbers, etc. ; however, ELD does not extract any 
data from TAB generated tables. Accordingly, if TAB is subsequently 
rc-executed to alter values of section properties, joint position coordinates, 
etc. , it is not necessary to re-execute ELD, provided the connected joint 
numbers, table entry pointers, etc. , of the elements are unchanged. 
Alphanumeric names used to identify specific types of elements used to model 
frame/shell structures are summarized in Table ELD-1. The formulations on 
which these elements are based are discussed in Volume II of the SPAR 
Reference Manual, 

Table ELD-1 ; Summary of One and Two-dim ensional Elastic Elements 

Name 
E21 

E22 

J23 


Description 

General beam elements, such as angles, wide 
flanges, tees, zees, tubes, etc. 

Beams of finite length for which the 6x6 
intrinsic stiffness matrix is directly specified. 

Bar elements having only axial stiffness. 
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Name Description 

E24 Plane beam. 

E25 Zero-length element used to elastically connect 

two coincident joints. 

E31 Triangular membrane, flat, aeolotropic, using 

TM constant-stress formulation. 


E32 

E33 

E41 

E42 

E43 

E44 


Triangular bending element, flat aeolotropic. 

Triangular membrane + bending element, flat, 
aeolotropic, using TM and TPB7 formulations. 

Quadrilateral membrane, aeolotropic, 
using QMB5 hybrid formulation. 

Quadrilateral bending element, aeolotropic, 
using hybrid formulation. 

Quadrilateral membrane + bending element, aeolotropic, 
using hybrid formulation. 

Quadrilateral shear panel, using hybrid formulation. 


For purposes of explaining card input, we will consider ELD to be 
comprised of an array of sub-processors, one for each type of element. 
The function of each sub-processor is to read input cards defining all of 
the elements of a particular type. A sub-processor is activated by an 
input card containing the name of the element type . Sub-processors may 
be called in any order. A typical input stream is shown below. 

@XQT ELD 
E43 

' ' 

{Data cards defining all of the E43 elements in the structure) 

E21 

(Data cards defining all of the E21 elements in the structure) 


ELD 


Execution of each sub-processor results in production of data sets having 
the following names (EXY represents the name of an element type, such 
as E21, E33, etc. ): 

DEF EXY - Basic (integer form) element definitions. 

GD EXY = Group directory. 

GTIT EXY = Group titles . 

When a sub-processor is executed, the resultant output data sets replace 
all data previously generated by the same sub-processor in any previous 
execution. * 

Elements of each type may be assigned by the analyst to separate 
groups . Within each group, each element has an identifying index number. 
The grouping of elements serves various purposes, which new users will 
learn with experience. It is almost always best to use many groups. 
Definition of the elements in group n is preceded by an input card in the 
following form; 

GROUP n* - - - title describing group n - - -. All elements 
defined by input cards following such a card belong to group n. Groups 
must appear sequentially (group 1, group 2, — ). If no GROUP card appears, 
all elements are assigned by default to group 1 . The order in which individual 
element definitions appear following a GROUP card determines the index 
numbers they are assigned within the group . 


* It is planned in a future release to provide an "update" mode of operation 
allowing deletions, additions, and corrections to be made. 
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Element definition cards have the forms shown below. 

All Z -node elements: 

Jl, J2 

All 3 -node elements: 

Jl. J2, J 3 

All 4-node elements: 

Jl, J2, J3 , J4 

It should be noted that the order in which Jl, J2 - - appear on the element 
definition cards is very important, since the orientation of element 
reference frames is defined on the basis of this information (see Pig. ELD-1). 
All element- related quantities, such as stiffness coefficients, stresses, etc. , 
are relative to these reference frames. 


Optional 


NETOPT, NET(l), NET{2) 


NETOPT, NET(l)i NET(2), 


NETOPT, NET(l), NET(2), 


In the above, Jl, J2, — indicate the joints interconnected by the first 
element defined by the card. If the optional network operation parameters, 
NETOPT, NET(l), NET(2), — , are given, they result in definition of 
additional elements, in accordance with specific rules which will be defined 
subsequently. 

The section properties, material properties, etc. associated with 
each element depend upon the values of various table reference pointers 
prevailing at the time the element definition card appears. These pointers 
are summarized below. 
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Table Default 

Pointer Value Associated Table 


NMAT 1 


Material constants, MATC. 


NSECT 1 

NOFF 0 

NNSW 0 

NREF 1 


Section properties. E44 elements 
refer to Table SB. All other 3 and 4- 
node elements refer to Table SA. 

E21, E22, E23, E24, and E25 refer 
to BA, BB, BC, BD, and BB, respectively. 

Beam rigid Unit offsets, BRL 
{2 -node elements, only). 

Nonstructural distributed weight, NSW. 

Beam reference frame orientation, 

MREF (2-node elements other than E25, 
only). For E25 elements, NREF points 
to an entry in tbe A LTRE F table. 

For 3 and 4-node elements NREF is not 
a table entry pointer. Instead, NREF is 
used to specify the direction of action oF 
positive pressure exerted on the element. 
If NREF=0, pressure exerts no force on the 
element. S£NREF=1, positive pressure 
acts in the direction of the 3-axis of the 
element reference frame (see Fig. FLD-1, 
below). If NREF=-1, pressure acts in the 
opposite direction. 



(1) All systems right-hand 
rectangular. 

(2 ) O rientation of 1 , 2 axes of 

2 node elements depends on the 
entry in MREF table indicated by 
NREF. 

(3) For 3, 4-node elements, -J3 lies 
in quadrant 1 of plane 1-2, 


3 



Figure ELD-1. Element Reference Frames 
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An example of the use of the table reference parameters is shown 


below. 



NMAT=4: NSECT=7: 29,13: 13,74: NSECT=8: 74, 16: NMAT=2: 74,92$ 


EU3 also accepts input commands which have the effect of modifying 
the data given on subsequent input cards. These commands are summarized 
below. 


Command 


Meaning of n ( n default is z 

MOD 

JOINT* 

n 

Add 

n 

to all joint numbers. 

MOD 

GROUP* 

n 

Add 

n 

to all group numbers 

MOD 

NSECT* 

n 

Add 

n 

to NSECT pointers. 

MOD 

nmat= 

n 


IT 

NMAT " 

MOD 

NNSW= 

n 


H 

NNSW " 

MOD 

NREF= 

n 


:t 

NREF " 

MOD 

NOFF= 

n 


II 

NOFF " 

INC 

NSECT* 

n 

See below. 

INC 

nmat* 

n 

it 

It 


INC 

NNSW* 

n 

II 

II 


INC 

NREF* 

n 

It 

It 


INC 

NOFF= 

n 

11 

ft 
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A frequently encountered situation is the one in which a substantial 
portion of a finite element model is comprised of a regular mesh, which 
can be automatically defined by a mesh generation command, but the 
section properties are so nonuniform that it becomes necessary or highly 
desirable for simplicity of input preparation to associate with each element 
a unique entry in the section property table. A common example is an air- 
craft fuselage, where cutouts, doublers, stiffeners, etc., result in extensive 
variations in section properties over a mesh that is topologically and geo- 
metrically regular. A similar situation involving distributed nonstructural 
weight also sometimes occurs; e.g. thermal protection material of varying 
thickness. The INC command is intended primarily to address situations 
of this type . 

The INC commands are used to cause the associated table reference 
pointers to be automatically incremented by n as successive elements are 
defined. This is especially useful in conjunction with network generation 
options, which are defined in detail on subsequent pages. For example, 
the element definition command 701, 702, 1, 50 causes fifty elements to be 
defined, connecting joint pairs (701,702), (702,703) - - - - (750,751). If 
section property table entries 101,102, - - -150 apply, successively, to 
these elements, the following input would suffice. 

NSECT= 101; INC NSECT=1: 701,702, 1,50$ 

All MOD and INC parameters are internally reset to zero after 

conclusion of each sub-processor execution . The MOD and INC commands 

are not cumulative. That is, MOD NSECT 7: - - - - MOD NSECT 3; - - is 

not^ equivalent to MOD NSECT 10, 
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Reset Controls 

Only one reset command is permitted: 

OUTUtB = Destination library for ELD output. The TAB-generated 
data sets must already be resident in this library. 

Core Requirements 

ELD rarely requires ore than 1500-3000 words of working core 
space, regardless of the size of the problem. 

Code Release Data 


Level 7 (UNIVAC, CDC), July 1974, coded by W. D. Whetstone. 


mjfr 


Two**«iode element fietVsarft. gcnei^ato re* 

:* A ' .* ■ — 1 ^ -. .--I ■■ — v yq Mi 1 1 t 1 — 1 " "■ ' ■ " 


I^ISIETOPTal , 

- ,V ? V- 1 

NI= <NET«it) (defaults 1J 5 

NJ= NET(2J ^de fault= 

JINC= NEY(3). 

fen plied seqnen.se: 

N1=J1 

IEgFF*JS~?l 
D© 200 J=1,NJ 
DO £00 Jsl, NI 
N2=NHIDIFF 

Define element connecting nod© Nl to NF2 
100 N1=N3 
200 N1=J1 4 J * JEMC 


Example: 

J1 

if 2 

Netopfc 

m 

m 

JSNQ 

9 

12 

a 

4 

2 

ISO 

9 (1)*’ 

12 

(2) lg , v 

,m , * 8 

(4) 

21 

169 

U2 

11© 

118 


121 

Ov , , <5J -- 



c? 

1 

1 

“nsr 



--■ " — ■ - — ■ " 

*The order in which the elements ora defined is indicated by the faum.be e 
enclosed in parentheses. The index Slumber identifying elements within 
each group are determined by the order? in 'which the elements are defined* 
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Two-node clenri $nfc network generators Continued). 


If NETOPT=2 , 

NI= NET(l) (default= 1 }, 

NJ= NET(2) (defaults Uv 

J1NC= NET(ii» 

Implied sequence: 

NU J1 
ID1FF= J2-J 
DO 200 J=1,NJ 
DO 100 1=1, m 

N2 = Jjl + XDIFF (except for closing element, when I=NI) 
Define element connecting N1|N2, 

100 NUN2 

200 NUJ1 + JAJINC 


Example: 







J1 

J2 


Nefcopt 

m 

NJ 

jINC 

9 

12 


2 

; r- 

4 

2 

100 


9 




109 


(4) 




(8) 


< 5 ’ 

IS 



. 12 118 



> 112 

\ 

s(3)^ 



\ 


4 


15 




115 
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Two-node element network generators (continued). 


If NETOPT=3 , 

NI= NET(l), 

UNC= NET(2), 

NJ= NET(3) (defaults ) 

JINC= NET(4). 

Implied sequence: 

N1=J1 

N2=J2 

DO 200 J=l, NJ 
DO 100 1=1, NI 

Define element connecting Nl to N2 
100 N2 = N2 + IINC 

Nl = J1 + JINC*J 
200 N2 = J2 + JINC*J 


Example: 


J1 

J2 

Netopt 

NI 

IINC 

NJ 

JINC 

4 

37 

3 

4 

10 

2 

100 
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Three -node element network cenerators 


If Netopt= 1 , 


NI= NET(I), 

NJ= NET(2) (default= 1). 

Implied sequence: 

IINC = J2-J1 

JINC= J3-J2 
Nl- J1 

N2= J2 

N3= J 3 

N4= J3-IINC 
DO 200 J = l, NJ 

DO 100 1=1, NI 

Define element connecting Nl, N2, N3 
Define element connecting N3, N4, Nl 
Nl= Nl + I1NC 

N2= N2 + IINC 

N3= N3 + IINC 

100 N4= N4 + IINC 

Nl= J1 + JINC 

N2= J 2 + JINC 

N3= J3 + JINC 

200 N4= N3 - IINC 


J2 J 3 


Netopt NI 


(7) / (9) / (11) 

103 / 104 X 
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Three-node element network generators (continued). 



It Netopt=2„ 


NI= NET(l) (NX must be greater than 1), 

NJ= NET<2) (default=l), 

JINC= NET (3). 


Implied sequence: 


IINC= J3-J2 
Nl= J1 
N2= J2 

DO 201) J=l, NJ 
DO 100 1-1, NI 
N3“ N2 + nNC 

Define element connecting Nl,N2,N3 
100 N2=N3 


Nl= J1 + J*JINC 


200 N2= J2 + J* JINC 


Example: 

Jl J2 J3 Netopt NI NJ JINC 

2 5 7 2 6 2 30 
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Three-node element network generators (continued). 


If Netopt=3, 

NI= NET{1), 

NJ= NET(E) (default= 1), 

JINC= NET(3). 

Implied sequence: 

HNC= J3-J2 
Nl= JI 
N2= J2 

DO 200 J«l, NJ 

DO 100 1=1, NI 

N3= N2 + HNC (except closure when I=NI) 

Define element connecting N1,N2,N3 

100 N2=N3 

NI = Jl + J*JINC 
200 N2 = J2 + J*JINC 


Example: 

Jl 

2 


J2 

5 


J3 

7 


Netopt 

3 


NI 

6 


NJ 

2 


JINC 

30 
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Four -node element network generators. 


If Netopt=l , 

NI= NET(l), 

NJ= NET(2) (de£ault= 1), 

NK= NET(3) (defaults 1), 

KINC= NET(4). 

jL 

Implied sequence * : 


IINC= J2-J1 
JINC= J4-J1 
DO 300 K=l, NK 
Nl a J1 

DO 200 J = l, NJ 
DO 100 1=1, NI 

N2= N1+I1NC 
N3= N2+JING 
N4= Nl+JINC 

Define element connecting N1,N2, N3, N4 
100 Nl = Nl+HNC 
200 Nl = J1+J*JINC 
300 J1 = Jl+KINC 

Example” : 


J1 

J2 

J3 

J4 

NETOFT 

NI 

NJ 

NK 

KING 

2 

3 

23 

22 

1 

2 

3 

2 

200 

62 

63 


64 

262 

263 

264 



42 

22 


(5) 

43 

(6) 

44 

(3) 

23 

(4) 

24 

(1) 

(2) 

2 


3 

4 


242 

222 


< U > 243 

(12) 

223 

(10) 

(7) 

(8> 


244 

224 


202 


203 


204 


Note: J3 must be present, although not used. 
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Four -node element network generators (continued). 


If Netopt=2 * 

NI= NET(l), 

NJ= NET{2) (default=l). 

Implied sequence: 

IINC= J2-J1 
JINC= J4-J1 
Nl = J1 
DO 200 J=1,NJ 
DO 100 1=1* NI 

N2 = Nl + IINC 

N3 = N2 + JINC (except closure when i-NI) 
N4 = Nl + JINC (except closure when I=NI) 
Define element connecting N1,N2,N3,N4 
100 Nl = Nl + IINC 
200 Nl = J1 + J*JINC 


Example: 


J1 J2 J3 J4 Netopt NI NJ 
1 11 12 2 2 6 3 


3 13 23 33 
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The SPAK thermal element library consists of conducting elements 


K21, 2-node 
K31, 3-node, 
K41, 4-nnde, 
K61, fi-nodo, 
K81, 8-node, 


2D 


convecting elements 

C21, 2-node 
C31, 3-node, 2D 
C41, 4-node, 2D 


and radiating elements 

R21, 2— node 
R31, 3-node, 2D 
R41, 4-node, 2D. 

Element definition data is input via TELD. 

node numbering and local reference frames are pictured below. 


Element 


x 



Bar elements can have zero lengths to 
simulate point sources, convectors, and 
radiators. Zero-length bar elements are 
defined by prescribing the same node point 
at each end or by identifying zero-length 
bars in processor TGEO. 


Triangular elements must be numbered as 
shown. 


TA1-5 

A- 31 


1 


2 


y 




2 


page is 

OEWOR QBAL1TV. 


Quadr ilat i ml elements must be numbered as 
shown ami slu aid be nearly flat. 


First two nodes set + X direction axes 3rd 
node sets +Y direction. 



i 


z 6 



Pentahedral elements must be numbered as shown. 

They should have faces that are very nearly flat 
(all points in a plane) . Nodes 4~b muse lie above 
the x-y plane in + Z direction. 1 




z 



Hexahedral elements must be numbered as shown. 
They should have faces that are very nearly flat. 
Nodes 5-8 must lie above the x-y plane. 


A 


Element shape and size parameters are adjusted via RESET controls in 
processor TGEO. 
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Table TAl- r i < -n.. -tit Dol ini turn U.il.i Si*t_s 


In the* following t.ibut.it ion, ngroup identifies a group of elements of a 
particular typo. The parameters nmnt, nfilm, and nrad point to the applicable^ 
tables of CONDuc t i on , CONVect RADI at ion properties. The j *s are 

connected points. 

q3? PO0®- 

11 1 e number of elements in a group is indicated by the TOC item NJ. Each 
data line corresponds to rhe'j~tli element within its group (j=l,..., NJ) . 


Data Set Name 


Hi 

N2 

JO 

NI* 

Contents of 

data 

line 



TED 

K21 

ngroup 

7 

nmat 

h 


area 



TED 

K31 

ngroup 

13 

nntat 


3 2 

j 3 

thickness 


TED 

K41 

ngroup 

16 

nmat 


V 

j 3 

thickness 

TED 

K61 

ngroup 

27 

nmat 

j i 

j 2 


j 4 j 5 

^6 

TED 

K81 

ngroup 

35 

nmat 

j i 

J 2 

j 3 

i 5 

j 6 

TED 

C21 

ngroup 

7 

nfilm^ 

tt 

j i 

j 2 

circumference 


TED 

C31 

ngroup 

13 

nfilm 

Jl 

j i 

^ 2 

j 3 



TED 

C41 

ngroup 

16 

nfilm 


j 2 

j 3 

j 4 


TED 

R21 

ngroup 

7 

nrad 


h 

circumference 


TED 

R31 

ngroup 

13 

nrad 

h 

h 




TED 

R41 

ngroup 

16 

nrad 


^ 2 

J 3 

^4 


All 

group 

titles for 

element 

type K21 

(or 

K31, 

etc.) 

will reside in 


'8 


named GT K21, produced by AUS/ALPHA. The first block corresponds to group 1, 
the second to group 2, etc. 


* The NI values are greater than the number of quantities actually input to 
allow certain system computed information to reside in these data sets. 

^ If nfilm is negative, its absolute value points to an entry in a table of FILM 
transfer COEFf icients. 
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Temperature dependent material pinperry tables are constructed with 
AUS/TABLE as shown below. Each data line defines properties corresponding to 
temperature T. Properties at intermediate times are determined by linear inter- 
polation. The number of temperature points (parameter NJ) may be different for 
each table. If only one data line is p.iven (NJ=1), properties are constant. 


Data Set Name 

— PM 

page is 
quality 





Nl 

N2 

N3 

NT 

Data Line Contents 


CONDuction 

PROParci.es 

nmat 

9 

T, 

P* c, an a 2 2 “33 “12 

“23 “ 3 i(input all 

CONVection 

PROPerties 

nfi lm 

2 

T, 

h 

values » no 
default) 

RADIation 

PROPerties 

nrad 

2 

T, 

£ 


Film transfer coefficients which are functions of geometry (e.g. variation 
along an aerodynamic boundary layer) can be defined in one data set. The TOC 
parameters are Nl=l and NJ=number of film transfer coefficients. 

N1 N2 Data Line Contents 

FILM COEFficients hj, h^ h^j 

The SI units for the above material properties are listed below. If a differ- 
ent homogeneous set of units is used, the Stef an-Boltzmann constant must be re- 
defined via an SSTA or TRTA reset control (Sections TA1.5, TA1.6). 

T * temperature, degrees Kelvin (K) 
p = mass density (kg/m 3 ) 
c s» specific heat (J/kg K) 

= conductivities (W/mK) , defined with respect to the element 

coordinates shown in Section TAI.l. Any conductivity not specified 
is assumed to be zero, 
h = film transfer coefficient (W/m 2 K) • 

£ = emissivity (dimensionless) 
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TELD FOR THERMAL ANALYSIS 


RESET CONTROLS 

def aul t 

NOTED 0 Unit on which to save TED el name ngroup 

data sets. If NUTED is not reset, no 
data sets are created 

LRTED 896 Record length of TED el name ngroup data sets. 


Thermal elements are defined in exactly the same manner as structural 
and fluid elements are defined. 

The table pointers NSECT. are used. for thermal definition. 


TABLE POINTER: 


default 

NMAT 1 Material properties pointer 

nmat for conduction properties 
nfilm for convection properties 
nrad fop radiation properties 
(See Table TAI.2, page TA1-8) 


NSECT 1 Entry in the appropriate table 

constructed via AUS/TABLE 
K AREAS (K21 elements) 

K THICknesses (K31 , K41 elements) 
C GIRCumference (C21 element) 

R GIRCumference (R21 element) 

These data sets must exist 
before ELD is executed. For 
example; if K21 elements are 
to be constructed, the data 
set K AREA must exist. 


$ If nfilm is negative, its absolute value points, to an entry in a 
tgble of FILM transfer COEFficients 
{See page TAI-8) 


5*1.3 Date -Set Constrijptttrs 


Several kinds of dat&tfets may S>S go exited or modified % the TAB&E, SYSVEG, 

r 

ELDATA ? and ALPHA sub$roces£f©rs. Examples are summarized in the following 
tabulation: 


Subprobassor 

Kind of Data Set Produced or Operated On 

TABLE 

TABLE.** Examples: Nodal pressures, tempera** 
tures c and TAB-produced tables. 

SYSVEC 

SYSVEC.* Examples: System vectors of applied 
forces^ motions,, and eigenvectors* 

ELDATA 

ELDATA. * Examples: Element-applied load data*© „ ft* 
pressures',, temperatures, and dislocations* 

ALPHA 

ALPHA.* Examples: Arrays of alpha character 
strings -used as load ease titles,, or to describe 
eigenvectors. 


In the following, XICXX is pny of the subprOoesscrr names : TABLE* SYSVE0* ELDATA 

or ALPHA. *3*0 initiate production of a nev? date s©§ named SI N2 jn3 a4 } the 
following language sequence te used: 


XXXJg|®ptionsd parameters): N1 ND n3 n4? «-=■ data sards «*■«■ 


To Initiate modification (updating) of aa existing date set-, the following language ils 
used: 


XXXX, tl( optional parpKfiet^J; 351 H? cS ji 4? $ata garde. »«** 


♦See Section 2.5. 
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Si QUALITY 

The output data set, K1 N2 n3 n4, is resident in the current data destina- 
tion library defined by the last OUTLTB Statement . If omitted, N2, n3, and n4 
assume default values of AUS , 1 , dnd 1 . 

The following characteristics of the Oper, U (update) mode of operation are 
noted: 

• The data set is actually overwritten in mass data storage, rather than 
modified and stored in a new area. 

• Existing data sets must not be extended (i.e., additional blocks cannot be 
appended in this mode of operation). 

The following example input indicates a typical application of a data set constructor. 

In this case SYSVEC is used to create a data set, named APPL FORC 88, comprised 
of two blocks. The first block defines load case 1 of static load set 88; and the second 
block defines case 2 of set 88. 

(SXQT AUS 
SYSVEC$ 

APPLIED FORCES 88$ 

3$ Direction 2 and 3 forces. 

List of joint numbers. 

*7.3$ 3fg, f«j at joint 7. 

8.3$ »' " 8. 

9.3$ t " " 9. 

CASE 2$ 

1- 1: J=3.0j 44.6$ fj at joint 10. O 

1^.3: J= 7: 19.2$ X3 at joint 7. 


CASE 1$ 

I— 2-, 

J= f, 9$ 

7.2 

8.2 
9.2 


Enter SYSVEC 
Name of output data set. 


5.T. 3-2 


Many examples of SYSVEC, etc. , input arr given in Section 6. It is suggested 
that these examples be scanned briefly before reading Sections 5. 1. 3. 1-4. 

Because of the central role played by the data set constructors, it is strongly 
recommended that new users perform a series of test executions of them, using 
DCU/PRINT and/or VPRT to verify that the intended results are achieved. 

The TRAN sub-activity of TABLE provides a very general method of trans- 
mitting information from one data set to another. New users should not attempt 
to use TRAN without advice from an experienced analyst. 


5. 1.3-3 
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AUO/ 

TABLE 

5. 1.3.1 TABLE. Execution begins as indicated below. If appropriate, any of the 
underlined words may be omitted. 

TABLE, UfNI- ni. NJ= nl) : N1 N2 n3 n4 : —data — 

The optional parameters NI and NJ have the same meaning as defined in 
Section 2.5; i.e., each block of the output data set is a rectangular matrix of dimen- 
sion (ni, nj) . The default value of ni is 1, and the default value of nj is the 
number of joints in the structure. The input defining the nth block of the output 
data set has the following structure: 

BLOCK n 

OPERATION = SUM $ or XSUM, or MULTIPLY, or DIVIDE 

J= jj; J 2 : igi - - - $ , a list of columns. * 

i 2 3 

ef ef e. $ Data record 1, applies to column j. 

J 1 J 1 3 1 3 1 1 

12 3 ni 

ef ef e. e. $ Data record 2 , applies to column j„ 

h h h h 2 

etc. 


*Loop-llmit format is also permitted (Bee Section 2. 3). 
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If there are fewer data records of the latter type than there are j’s in the 
list of columns, the last record will be associated with each of the remaining j’s . 
The moaning of the ej's depends upon current arithmetic mode, as established by 
the OPERATION 3 XXXX statement. In the following, is the element in row i 
and column j of block n of the output data set. The range of i is 1 , 2 , ~ “ ~ ni, 
and the range of j is Jj_ • Jo . 3^ • 


Operation 

Arithmetic mode 

SUM 


- *J + e \ 

XSUM 

■i 

-e* 

MULTIPLY 

*s 

~ x ) e J 

DIVIDE 

■5 

- x W 


The sequence (OPERATION-XXX List of j’s data records for 

each j in the list) may be repeated as many times as necessary to produce the desired 
data in block n. 

The OPERATION statement continues in effect , not only for the current block 
but for all subsequent blocks as well , until another OPERATION statement supersedes 
it. The default OPERATION is SUM. 
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The following control statement is frequently useful: 


1= i* i^ — i m $ (each i^ between 1 and ni) 


This statement indicates that the data records will have the following content: 


For example: 


j 1 / i m , 
6 J * 6 j — e J * * 


BLOCK n: etc. : - - - $ 


1= 3,5: J= 7,9$ 


3 5 3 5 3 5* 

e,.e 7 > V e a : VV 


1=1: J= 10: e 10 $ 


1 2 

The 1= i , 1 — statement continues in effect, not only for the current 

12 

block but for all subsequent blocks as well, until another 1= i , i — statement is 
encountered. Default is, in effect, 

1= 1,2, — ni , 
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A US/ 

TABLE 

The following statement may immediately precede a list of j f s: 

DDATA= dj.d,, -d m $ 

When data records are being duplicated, the DDATA command has the effect 
of incrementing each item in the data record by the corresponding dj . The DDATA 
statement will apply only to the j— list immediately following the statement. 


Example: 


1= 1, 3: DDATA= .2, .1; J= 2, 5: 7., 8.$ 
The foregoing statement produces the same result as: 


1= 1,3 : J= 
7 . 0 , 8.0 
7.2, 8.1 
7.4, 8.2 
7.6, 8.3 


2,5$ 


$ 

$ 


$ 



e 


3 

2 


e 3 * 6 3 


e 4’ e 4 


V e 5 
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To cause an integer constant, jshift, to be added to each J in subsequent 
lists of j'a, the following command is used: JSHIFT* jshift. The jshift parameter 
will remain in effect until superseded by another JSHIFT command. 

Rules for ordering control and input statements: 

• No statement may precede the first BLOCKS n statement 
(Exception: if operation is confined to block 1, the BLOCK- 1 
statement may be omitted). 

• When the BLOCK= n statement appears, block n - 1 must already 
exist. 

• No control card of any kind (OPERATION = — , 1= — JSHIFT* 8 — , 
or DDATA= — ) may appear anywhere between the beginning of a 
J-liflt and the end of the subsequent data records ^the ej's^ • 

Synonyms. If desired, the control word CASE may be used instead of BLOCK, 
and JOINT instead of J. 

Multiple data records on a single card . . In the preceding, it has been 

indicated that each data record should contain exactly ni words, or m words if the 
12 m 

statement 1*1,1 — i m is used. It is also permitted to place several such 
complete data records on the same card; i.e. , both of the following are legal cards 
and have the same effect: 

1.2, 2.9: 4.5, 2.7: 3.2, 9.$ 

1.2, 2.9, 4.5, 2.7, 3.2, 9.$ 

An error stop will occur if a card does not contain complete records. 

5. 1.3. 1-5 
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LID 


/» 

The TRANSFER Statement. This command provides a general method of 
extracting real data from designated areas of one data set, and adding them to data in 
designated areas of another data Bet. TRANSFER statements may appear as the first 
or last statement In a TABLE execution, or immediately prior to any BLOCK state- 
ment. The command form is as follows: 

TRANSFER( List of arguments P| ,P= Vg t — ) 

Any or all of the arguments summarized In the following table may be used. 


Argument Default 

Name Value Meaning 

SOURCE None Name (a single alpha word) of the source 

data set. May be established by a 
DEFINE statement. 

LI 1 Successive blocks of die source data set 

are added to blocks LI, LI + 1, — L 2 
of the destination data set unless Ll=&2. 

L2 1 

IUM Leaser of b B * See Fig. 5. 1.3. 1-1 

JLIM 1 

SSKIP 0 

DSKIP . 0 

SBASE 0 

DBASE 0 


The TABLE/operation statement has no meaning In TRANSFER. 
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AUS/ 

TABLE 


Source 
Data Set (S) 


Destination 
Data Set (D) 


T 


SBASE 


r 

■ J 

DBASE 

- 4 - 

„ \ 


TTJM 

i 

\ 

T 

DSKXP 

ILEM 

i — H 

J ( 

* 

I 

Block 

T 



Fortran equivalent, 
repeated for blocks 
LI, LX + 1....L2: 

IS -SBASE 
ID - DBASE 

DO 20 J - 1, JLDI 
DO 10 I - 1, HIM 
IS - IS 4* 1 
ID - ID 4* 1 

10 D(!D) - D(xD) + S<IS) 
IS-IS+ SSK3P 
20 ID - ID + DSKJP 



Fig. 5. 1.3. 1-1: TRANSFER Parameters 
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Figure 5. 1.3. 1-1 illustrates the function of the TRANSFER statement for the 
case in which the number of blocks in the source data set is equal to the number of 
blocks in the destination data set (i.e. , L2-L1+1). 

If Ll= L2, the function is similar, except that successive blocks of the source 
data set are transmitted to successive areas within block LI of the destination data set. 

The following are a few examples of the TRANSFER statement. 


&XQT AUS 

DEFINE OLDV= 2 VIBR MODE 200 1 1, 10 $ 

TABLE (NI=6, NJ=1700): VIBR MODE 11 $ 

TRANSFER ( SOURCE^ OLDV, 11=1, L2=10) $ 


Salvage vibrational modes 
from a previous analysis for 
use as initial approximations 
of modes of a modified struc- 
ture having additional joints. 


DEFINE U7= STAT DLSP 117 $ 

TABLE (NI=3, NJ=400): DEFORMED POSITION $ 

TRANSFER! SOURCE^ JLOC) $ 

TRAN(SOURCE*U7, ILIM=3, JLIM=400, SSKIP“3) $ 


Store node point locations of 
the deformed structure in 
DEFO POSI. 

JLOC= original position. 

U7= Joint motion due to static 
load case 7 of set 1. 


DEFINE X= VIBR MODE 12 18 $ Form a single-block matrix having 

8 columns, each column corres- 
TABLE( NI=30, NJ=8); T MODE 12 $ ponding to a mode. Each column 

contains all 6 motion components 

TRANSFER! SOURCE=X, SBASE=600, ILIM=30) $ of joints 101-105. 
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r» 2 DCU - DATA COJUPLKX UTILITY PROGRAM 

Function - DCU performs an array of utility functions, as subsequently 
defined in this section. New users should review Sections 2.2 and 2.5 before reading 
this section. 

In defining DCU command syntax, the symbol Lib will be used repeatedly 
to represent a SPAR library internal designation (1, 2, etc.), and Id will identify 
one or more data sets . Id may be in any of the three forms indicated below, unless 
specifically stated otherwise: 

(1) The four- word data set name, MASK- filled. 

Example: K SPAR is the same as 

K SPAR MASK MASK. 

(2) An integer, n, indicating the data set associated with 
sequence number n in the table of contents. 

(3) Integers n, m, indicating (m-n 1 ) data sets associated 
with sequence numbers n, n+ 1, - - - m. 

(4) Omitted - meaning all data sets in a library. 

Consider as an example, the most commonly executed command, 
TOC Lib Id. 

Examples: 

• TOC 1 K SFAR$ Id form (1) causes a single line of the 

Table of Contents of library 1 to be printed, i. e. the line 

5.2-1 
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net: 


corresponding to the first data set named K SPAR MASK 
MASK to be located. 

ft 

• TOC 1 27$ Id form (2) causes printout of line 27 of the 

Table of Contents of library 1 . 

* 

q TOC 1 32,50$ Id form (3) causes printout of lines 32 through 

50 of the Table of Contents of library 1. 

• TOC 1 $ Id form (4) causes the entire Table of Contents of 
library 1 to be printed. 

Other commands currently available in DCU are summarized below: 

9 DISABLE .Lib Id$ , Data set(s) are marked as disabled. 

The data set(s) are still present in the library, but cannot 
be accessed until they are re-enabled via the following command: 

4 

9 ENABLE Lib, Id$ . Only forms (2) or (3) jd £ Id are allowed 
for this command. 
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• PRINT Lih ld$ , or 

PRINT l.ili Nl. N2, n3, nl, jj, j >( , i^, i„. bj, b (> . The one or more 
identified data sets are printed in tabular form*. Only matrix-form data 
sets (e.g., TABLE, SYKVEC, ELBATA, and ALPHA forms, as defined 
in Section 2.5) should be PRINTed. If the second form is used, the 
printout wall be restricted to columns jj through jg, rows (items) 
i ^ through i 9 , for successive blocks b^, 1 — bg . Using the 

terminology of Section 2. 5, default parameters are = 1 , jg = NJ , 
i. - 1 , i = NI, b, ~ 1 , b 9 - total number of blocks in the data set. 

NOTE 

PRINT displays data set items sequentially, 
so that matrices appear in transposed 
form . 


Data produced by processors E , EKS s TOPO , K , M , KG , and INV 
should not normally be named in a PRINT command. 

• COPY Libj , Libg , Id$ . Copy the indicated dgta sets from Lib^ 

to Libg , Disabled data sets are not copied. This is the recommended 
method of packing libraries. 

• XCOPY Libj , n, Id$ . The indicated data set (in Lib) is written on 
ordinary sequential file n in a sequence of physical records identical 
to individual bio cits of the data set as it resides in Lib . As an example. 


* An auxiliary command, NCPL= n, controls the number of columns 
per printed line. Default NCPT.= 10. For teletype display, select 
NCPL- 5. 
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I- V. 


X®0P? %} 5, •anuses ibdgsowectojo ?b»&e «aritteij •iQiitj? 

a file known bxtefbaay as SPAR-E fllMStfAQ -*r B&A$m (CDC&b 
corresponding to 3i - 5» She output file will .ponfcaiJi /one ©hypipai 
record for <feach eigenvector, 


• XLQADn , Lib nwords . ni njnj , „. tVPe,.,„ .. Nl ,,„ _NSL . -i)#,, ,-n4 : $ ., 'fhiS 
con^naii^ causes data $som sequential Hie % to he (loaded .as ii'dafg./sefc 
named N2 , $3 , n4 in Lib , The e> their 5 parameters have the 
same meaning;^ defined in Section SJ? !9 As. jm ■©sample^, ?§upposq & 
sequential filc% SPAR-D^ contains five blocks {(physical £ece>rds 
produced by direct binary writes, f not r.nfoi»mut£ed Fortran Writes^ 
of real ilata and that each block is a matr&c with '6 row® and 100 columns. 
To load these datfi into library 1 as a data set named XX YY 1 
the following command would be used: 







°D 

°o 


4 


o 


o 

o 






ai 


o 

o 


XLOAD 4,1, 3000, 100^ 60&, XX YY 2. 2 $ 

• REWIND n $ . Used in conjunction with XOQPY x XLOAEJj- ithiS cam<* 

mand iiauses sequential file n to he rewound set to. starting 

point). It should be noted that neither IfcCOPY nor- XL.t)AD rewinds 
sequential files either befdrif) stftet Ih© ?date> transmits 10% $«jt that 
one sequential file may contain aomy multiblock data sets,, 

• SCALE Cg$. TMs command may precede a be' 

I upp 

If it does, the output data ##ts will = (so^cc^ >> ©g 4 
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• CHANGE Lib Id ,,, Id S. This command causes the 

. old new T 

name of data set to be changed from to ^ new * 

the full 4 -word name form is permitted for both Id's. 

fr 

• DUPLICATE Lib^, Lib^. Lib, is creatod identical to 
existing Libj, including disabled data sets, if any. 

• TWRITE Lib$ . Lib is written onto tape nt (see NT APE 
command). The complete library is written in physical 
records as large as the available core will allow,. 

9 TREAD Lib$ . Lib is read from tape nt (see NT APE command). 
Available working core space must be as large as it was when 
the T WRITE was executed. 

• NT APE)* nt$. The internal unit number of the tape to be 
used in the next T WRITE or TREAD command is nt 
(default* 20)., Note that logical 20, is Univac file SFAR-T, 

§DC file SPARLT. 

• STORE Lib, 1D$. Lib is stored as a dat^ set named Id, 

in library nZ (see LIB LIB command). Id may only be a full 
4»ward name, the first two words of Which are typeless. 
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$ RETRIEVE Lib, Id$ . The data set Id is recovered from the 
library njJ and constituted as library Lib. 

i 

• LIB LIB a n£$> The internal unit number of the library 
library is n£ (de£ault=12). Note that logical 12 is Univac 
file SPAR - L , CDC file SPARLL. 

* 

O TITLE Lib 1 - - - Alphanumeric title for Lib - - - 
The label-field title is embedded in Lib, and will be 
displayed at the beginning of each table of contents 
printout produced by a TOC command, 

• STATUS Lib 3 . The number of library entries and the 
current I/O counts for Lib are printed. 

® ABORT n$ . To cause an error -abort if an abnormal 
event occurs in DC U, setABORT=l$. 

Core Requirements. Working core must be sufficient to accommodate one block 
of each data set transmitted through core (e.g. via COPY, XCOPY, etc.). See 
also the discussion of TWRITE and TREAD. 

Code Release Data. Level 9, July 1975, coded by W. D. Whetstone. 
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(w/o updates for LaRC installed Jan. 1978) 
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ABSTRACT 


The GEOMPLT program is an interactive graphics program for finite 
element and thermal network models. The program was developed for use on 
a Tektronix 4000 series CRT thermal connected to a UNIVAC 1100 series 
or CDC 6600 computer. The program has the following features: 

1) Plot data frcm the NASTRAN, SPAR, STRUDL (UNIVAC only) and 
CINGEN progress. 

2) Automatically expands core for NASTRAN and STRUDL bulk data 
(UNIVAC only). 

3) Create hard copy plot . le for drum plotters. 

4) Create input data files ~>r CONTUR program. 

5) Provides an alternate print file for post plotting diagnostics. 
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editor routine. Mr. J. Chiou was responsible for the SPAR, conversion routine. 

Mr. W. E. Shultz and Mr. L. D. Craig implemented the thermal network and 
XY-Plot capability. 


* 
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1.0 INTRODUCTION 


1 

It has long been recognized that a graphic presentation of a finite element 
model is necessary to validate the model prior to structural analyses. Most finite 
element programs contain some form of general purpose plotting. The plotting 
capability of such programs as NASTRAN (Ref. 1) and other general purpose structures 
programs are quite extensive and very useful; however* they generally require large 
core and offline plotters which results in delayed processing and, therefore, ; 

delays in data editing. With the increased use of mesh generators (Ref. 2) operated 
in a demand terminal environment, it has been necessary to develop stand alone, 
small core, graphics programs for use in this new modeling environment. ; 

The GEOMPLT program in its current state of development, provides a program 
which operates with small core (as little as 22K), on a Tektronix CRT graphics j 

terminal connected to a UNI VAC 1100 series or CDC 6600 computer. The initial development 
of the program began in 1968 with a small off line plotting program which utilized a 
small core computer (16K) and a Cal Comp drum plotter. The program was first used 
in conjunction with the SAMIS program and since that time it has been modified s 

and improved. In 1970 it was converted to plot NASTRAN data and to operate on 
the UNIVAC 1108 computer system. In 1971 an interactive graphics version was 
written for use on the UNIVAC Advanced Graphics terminal. In 1972 a version for 
use on the 4000 series Tektronix CRT terminals was completed. Since that time, 
the program has been expanded to plot data from other finite element programs and \ 

conductance/ capacitance networks for thermal analysis programs. j 

1 

. J 

The program is written in FORTRAN with the exception of two assembly routines 
(UNIVAC only) used to control system parameters. The program in its present configuration j 
and combined with the use of mesh generators, provides the necessary software for 
complete model generation, and viewing. All data can be prepared and checked 
prior to execution from a demand terminal. 

The following sections of this manual describe the program installation and 
initiation, program features and their use, and finally a brief description of 
each subroutine. 


2.0 OPERATING ENVIRONMENT 


The most recent version of GEOMPLT program (Version 4*0) is designed to 
operate on the UNIVAC 1108 Exec 8 Level 31 system, with a Tektronix 4000 series 
CRT terminal. Earlier versions exist for the CDC- 6600/Tektronix 4000 series 
terminal, UNIVAC 1108/CalComp model 565 drum plotter, UNIVAC 1108/UNIVAC Advanced 
Graphics terminal and IBM 1130/CalComp Model 565 drum plotter. 

The earlier versions represent historic developments of the program leading 
to the current version described in this manual. 


2.1 PROGRAM INSTALLATION 
UNIVAC PROGRAM SETUP 

The program is contained on a 9 track tape, recorded on three files in 
800 bits per inch (BPI) density. A fourth file contains a small NASTRAN bulk 
data deck for testing the program. 


All files were written with @C0PY,GM executive instruction. Therefore, 


they may be read into the local system with the .following instructions: 


©ASG,T IN,T,# 

@ AS G, UP ABS,F2 
• 0ASG.UP SR, F2 
@ASG»UP TCS,F2 
©ASG.UP 4,F2 
©REWIND IN. 
8C0PY,G IN. ,ABS. 
0COPY.G IN., SR. 
0COPY.G IN. ,TCS. 
0COPY „G IN., 4. 
©FREE ABS. 

©FREE SR. 

©FREE TCS. 

©FREE 4. 

©FREE IN. 

©FIN 


• GEOMPLT PROGRAM TAPE 


(Absolute element => .MPA) 

(Source and Relocatable element) 
(Terminal Control System Library) 
(NASTR/N Bulk data pimple) 
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2.2 PROGRAM EXECUTION FROM TEKTRONIX TERMINAL 


The GEOMPLT program as described In this manual is currently operated by 
Sperry in two environments. 

The first uses a Tektronix 4013 connected to a TTiJIVAC 1108 

3x2 system, Exec 8 Level 31, with 1200 baud transmission rates* On this system 
the program is cataloged using the SECURE processor. The user enters the following 
commands to initiate a plotting session: 

Al« Set modem to "talk". 

A2. Dial (241), when high pitched tone is heard, push "data" button and 
set phone into receiver. 

A3. Clear screen - push "reset” 

A4« Type site ID - |pTK230| 

A5. When system responds, type run card j@RUN - - - -[ 

A6. Set qualifier for secured file iQQUAL SPERRY) 

A7. Secure program 
@SRG 

SPERRY*GE0M,* 

OUSE TPE$.,*GE0M. 

A8. Data to be plotted is obtained from one of the following sources: 

• A previously written tape 

• A previously cataloged data file 

• A mesh generator program 

• The C1NGEN program 

' • The TOPOLOGY I program 

• The SPAR program 

• A newly created data file using the DATA,ELT or ED processors and 
keyboard input. 

Data must be placed in or attached to file 4 for reading by GEOMPLT, 
e.g., @USE 4., BULK. (Data on file BULK), 
or 

@ASG,A 4. 


(Data cataloged o:i file 4) 



An exception occurs if the data to be plotted is from the SPAR program. 
In this case the TAB and ELD processors are executed and the print file is 
BRAKPT to a cataloged file 9. Pile 9 is the GEOMPLT input file for plots of 
SPAR data, e.g.,: 


0ASG.UP 9,F2 
OBBKPT PRlNTS$/9 

@XQT TAP 

' 

(SPAR INPUT DATA) 

a 

@XQT ELD 

a 
a 

(SPAR INPUT DATA) 

6BRKPT PRINT$ 

©FREE 9. 

@ASG,A 9. 

el n i n . i ■ — > 

A9. Execute GEOMPLT Program 


@ASG,T 

21,F2///500 

@XQT 

.MPA 


A10. The program is conversational from this point on and the user must 
respond to requests by the program in a step-by-step manner. 

The second system is a UNIVAC 1108 1x0 

system. Exec 8 Level 31, with a 300 baud transmission rate. On this system the 
program is stored on tape and must be cataloged by a batch job prior to initating 
the graphics sessior 

@RUN 

@ASG,T IN,T, delivery tape jj 

@ASG,UP GE0M.F2 

@ REWING IN. 

@C0PY,G IN. ,GEOM. 

0FREE GEOM. 

@FIN 
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( 


With the program cataloged, the Graphics Job can then be started: 


Bl. Turn Acoustic Coupler modem to ”on" 

B2. Dial UNIVAC phone number, when a high pitch tone is heard, place phone 
into Acoustic Coupler. 

B3. Type Site ID (switch on 4013 to "LINE") 

B4. When the system responds, type in a run card, 



B5. The requirement for data files 4 or 9 are the same as described in 
steps A. 8 or A. 9. 

B6. To execute GEOMPLT program 

@ASG,T 21,F2///200 
OXQT GE0M.MPA 

B7. The program is conversational from this point on and the user must 
respond to requests by the program in a step-by-step manner. 


f. 


i. 
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CDC PROGRAM SETUP 


A. Computer Configuration 


1. 

Computer 

CDC 6600 

2. 

Core size required 

26K Decimal 

3. 

Language 

FORTRAN IV 

4. 

System 

SCOPE 3.4 or NOS 

5. 

Plotter required 

Yes - Tektronix model 4000 series 

6. 

Card punch required 

Yes - If punch option on 

7. 

Tape Assignments 

None 


B. Estimated Running Time 

CPU time in GEOMPLT was 7.03 seconds with a problem size of 221 elements 
and 411 grid points, and 1 plot. An additional 86 seconds was required 
when data was generated for the CONTUR program (CONTUR =1 or 9). 

C. Restart Procedure 

Not required due to small amount of run time required. 

D. Deck Sequence 

The finite element data is read into the CDC remote system terminal and then 
batched to Tektronix graphics terminal. 

For Intercom Job: 

Notes : 1 . A=send button 

2. Comments are enclosed in brackets { } 

3. SS = two numbers assigned by the operating system 

4. v=Blank space 

UT-2QU Terminal (LX) 

R, BULKa 

BATCH, BULK, PRINT, LW, ZZa 

Tektronix Terminal (LW) 

H, 0a 
EA 

BATCH, IZZLXSS, LOCALa 
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BATCH, IZZLXSS, RENAME, TAPE4 

REWIND, TAPE4 

ETL = 500A 

ATTACH, NASB, LPXXXA2, ID = LPXXX, CY = 2a 

ATTACH, AGII, LPTEKTRONIX , ID = LPPATRICK, CY = 2A 

XEQA 

LDSET, LIB ® AGIlA 

NASBA 

TITLEa s Up to 80 characters may be used for the title] 

v $JNPllTv lEnter options desired from Table II, for example, ENO = 1, 
GNO = 1 , PV = 1 , etc}$A 

v $PVIEWv{Use only if PV = 1, in $INPUT, for example, 

JPV(l) = first element number, - 0, last element number 
of a desired set of elements to be plotted, first element 
number of another set to be plotted at the same time, -0, 
Last element number of the second set of elements to be 
plotted, single element number, single element number, 
etc}$A 

EDITA (Option after each plot if editing of NASTRAN connection 
or GRID cards are desired.} 


or 

Na (No editing at this time and the next plot is initiated 
by typing title, $INPUT and $PVIEW.} 


45 (5-10-74) 


B-10 


3.0 PROGRAM FEATURES 


The GEOMPLT program performs the following functions! 

• Graphically displays finite element models (NASTRAN, SPAR or STKUDL), 
including partial views, rotated views, a resistance network 

with dashed lines and element and grid point numbers. 

• Creates data files of geometry and connection data for use by the 
CONTUR program (Ref. 3). 

• Edits finite element data through addition, deletion, and modification 
of bulk data card images. 

• Creates an alternate print file which provides a permenant record 
of the terminal session for review at a later time* 

• Writes a data file for obtaining hard copies of plots via CalComp 
or Houston Instrument drum plotters. 

• Automatically expands core to the minimum required for thB NASTRAN 
problem. 


3 . I GEOMETRIC PLOT 


Purpose : To display a 3-D view of the fi»t$e rapdeS a $erra|ha',i £KT? 

Options : • Plot NASTRAN, C INGEN, SPAR ®r STRUDL (UMlVAC ontyfc model’s., 

• Print grid* element, Capacitance and conductance numbers., 

• Draw any element type listed in Table I* 

• Plot partial views to obtain a @iear view of -a®c&3 areas £j# sejeq£er# 

elements. 

• Print a title at the bottom of the frara®. 

• Change plot scale factors for increasing, ©r flecx<sasing the size 
of the image plotted. 

• Rotate the view about three ones. 

• Interchange axes to obtain a rotated view or view utilizing a 
left handed coordinate system. 

Procedure : The various plot options are selected via. the Namelist's $ INPUT 

and $FVIEW. Table II contains a list of the $INEUT Name list parameters,,, 
a description cf the parameters, and the default value set ’by the 
program at the start of the program. All of these parameter's may be 
omitted from the Namelist Inpat or they may be <?t»angec£> prior to each 
plot. Parameters marked by * in Table II retain am Input value until 
• changed by a new input value; ail other parameters reverb to their 
default values once the plot is complete. 

The request for Namelist parameters is made conversationally by the 
program, l.e., a message appears on the screen requesting that the 
appropriate Namelist be supplied by the terminal operator at the 
appropriate time. 


The partial view is selected via the $PVIEW Namelist. Element numbers 
are listed thru strings of numbers and/or the use of (-0) symbol td 
indicate consecutive numbers, e.g. , $FVIEW JFV(l) Q 1,3,5,-0,10,13 $ 
is the same as, $PVIEW JPV(l) » 1,3,5,6,7,8,9,10,13$ 
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Note - Namelist Input roust start In Column 2 and end with a 


Parameters may be listed on several lines by using a 'V at the end 
of one line and beginning the next line in column 2 or greater. 

To terminate the program type "N" after end Of current plot. 

limitations : • Only elements listed in Table X may be plotted. 

■ The number of elements plotted depends upon the size of core 
available. Since the program will dynamically expand core the 
~ size limit is determined only by local restrictions on using a 

Tektronix terminal. 1305 grid points and 1305 element with an 
average of 4 connections per element will fit In 36K core. 

* 100 plots are permitted as a default.. If more are desired set 

a 

Namelist parameter NPLOT>IOO. 

Core Size Calculation : 

CORE « (IBANK + DRANK) + OPEN CORE 
1BANK+DBANK - 21651 

OPEN CORE « 4 x NUMGRD + 4 x NELEM + JSIZE + 10 + NJPV 
NUMGRD - Number of Grid Points 
NELEM « Number of Elements 

JSIZE ■ Number of Connections ■ E (NELEM x Number of Grid Points 

per Elements). 

NJPV - (300-NELEM/2) , if NELEM -600 
* 0 , if NELEM >600 
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Table I. List o£ Bulk Data Cards Used by GEOMPLT tj 


CBAR 

CTRIA1 

CONROD 

CTR1A2 

CROD 

CTRMEM 

CTUBE 

CTRPLT 

CTORDRG 

CTR1ARG 


Element Types Blotte d 
CQDMEM 
CQDPLT 
COUAD1 
CQUAD2 
CSHEAR 
CHRIST 
CTRAPRG 


CTETRA 

PLQTEL 

C WEDGE 

CHBDY: 

CHEXA1 

POINT 

CHEXA2 

ALINE 

C1S3D8 

REV 

C1S3D20 

AREA3 

AREA4 


Other Types Used 


GRDSET 

GRIB 

CORD2R 

CORD2C 

CORD2S 


All other types ignored 


Table 11. Namelist $INPUT Parameter Definition (UNIVAC only) 


Parameter 

Description 

Default 

Other 

baud* 

Transmission rate from computer to 
terminal, BAUD/ 10 

120 

Integer 

CONDNO 

Print conductor numbers for thermal 
network 

0-N0 

1-YES 

CONTUR 

Saves XY plotter coordinates on file KC 

O-NO 

1 or 9-YES 


and connector data on file KD to be read 
by CONTUR program. NOTE: When CONTUR = 1, 


the following default changes are used 
KTEST = 1 
KD " 11 
HP.IOT = 0 
NFLOT = 1 


ENO 

Print element numbers at element centroids 

0-N0 

1-YES 

GNO 

Print grid numbers at grid point 

O-NO 

1-YES 

IHUSTN 

Write hard copy ^lot commands on file 25; 
request to permanently save plot on file 
26 is made afte*. plot is complete 

O-NO 

1-YES 

KC 

File number to use when CONTUR = 1 

10 

Integer 

KD ’ 

File number to use when LTEST or KTEST = 1 

7 

Integer 

KTEST 

Option for CONTUR element punch option 
if KTEST = 9; plotting is supressed 

O-NO 

1 or 9-YES 

LTEST 

Option for CONTUR XY data punch options 
if LTEST = 9; plotting is supressed 

O-NO 

1 or 9-YES 

HBD0T 

Option for multiple plots. Set MBLOT ** 0 

1-YES 

O-NO 


to terminate plotting 



NCHBDY 

NASTRAN thermal elements to be plotted 

O-NO 

1-YES 


NODENO 

Print capacitance node number for thermal 
network 

O-NO 

1-YES 

a 

NPLOT 

ftimber of plots if MPLOT - 1 

10 

Integer 

i 

KTEST 

Option to print rotated coordinate data 
on file 21 

O-NO 

1-YES 


c 
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Tabla 11. (Concluded) 


Parameter 

Description 

Default 

Other 

PROGRH 

Type of input data on file 4„ 9 or 30 

1 = NASTRAN (file 4) 

2 * C INGEN DATA only plotted (file 30) 

3 » CINGEN and NAS IRAN element data 

1 

1-5 


plotted (file 4 and 30) 




5 ” Snm (HtaV ( UNIVAC 



PV 

Option for partial views; if PV - 1 
Namelist $PVXEvJ is used to specify 
element numbers to be plotted 

0-N0 

1-YES 

SCALEX* 

Option for adjusting the apparent width 
of the screen 

10 

REAL 

SCALE?* 

Option for adjusting the apparent height 
of the screen 

7 

REAL 

STUB 

Option to place a symbol at the grid point 
29 “ • 

29 

Integer 

X* 

Angle of rotation about x-axis (deg) 

34.27 

REAL 

Y* 

Angie of rotation about y-axis (deg) 

23.17 

REAL 

2* 

Ang’e of rotation about z-axis (deg) 

43.0 

REAL 

XAX1S 

Axis to be placed in screen x-axis 
C-f left to right) 

1 

1-6 

YAXIS 

Axis to be placed in screen y-axis 

2 

1-6 

9 

<+ up) 



ZAX1S 

Axis to be placed in screen s-axis 
(+ out of screen) 

3 

1-6 


Structure yes may be interchanged with screen 



axes for changing viewj, where: 1 =» x 9 2 = y 9 
3 B s 0 4 = -it, 5 ~ °y 0 6 “ -s 




NOTE: Parameters marked by * retain input value until changed by a new input 

value. 
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3.2 DATA EDITOR 
( 

The data editor present in Level 3.0 has been removed in Level 4.0. This 
was done to minimize the core requirements of the program and thereby Increase 
the problem size for a given core size. Those who wish to edit their data may 
do so using the Standard Text Editor available on the UNIVAC system. When 
editing is complete, the revised data may be viewed by re-executing the GEOMPLT 
program. 



— f\ 
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3.3 HARD COPY 


Purpose : To obtain a high quality hard copy of the plot displayed on the screen. 

The "hard copy" is drawn on a drum plotter of the CalComp or Huston 
Instruments type. 

Options : The user must select the hard copy option prior to plotting a view, then 
decide to keep the plot or discard it after the plot is complete. 

Procedure : The user sets parameter IHUSTM » l at the time NAMELIST $ INPUT is 
specified. 

« The program then writes general purpose plot commands on file 25 
While the plot is being generated on the screen. 

■ When the plot is complete the user is given the choice of savlr.g 
or discarding the plot. If he elects to save the plot, file 75 
(scratch file) is automatically copied to file 26 (save file' and 
file 25 is rewound for the neizt frame. 

• The file created by GEOMPLT will not drive a drum plotter directly. 

A special program must be written to read file ^6 and write ttu. 
proper plot commands, depending upon the type of drum plotter avail- 
able at the local computer site. File 26 contains the following 
information : 

First record: Title Format (13A6„A2) 

Second record: Krain£mum 9 Yminimura, Xmaximum, Ymaximum Format 

, (4E15.7) 

Following records: BRANCH, X a Y 9 ID Format (I5,2E15.7,I6) 

Where branch indicates the appropriate plotting action to be taken: 

1 = Not used 

2 => Draw a line from present pen position to the X,Y 

coordinates 

3 « Move pen to the 3E 0 Y coordinates 

= Write the Grid ID at the X,Y coordinates 

5 « Write the Element ID at the X,Y coordinates 

6 a Plot is finished, read a new first record 


7 => A symbol should be written at the X,Y coordinates 

8 » Write an X at the X,Y coordinates 

9 « Write a Y at the X,Y coordinates 
10 = Write a Z at the X,Y coordinates 

Last record: Title Format (13A6,A2) where the first word in the title 

is ENDRUN 
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3.4 COHTUR PlJOT 


Purpose ; To create two data files UO & H) which contain transformed grid 
coordinate data and connection data for use by the CONTUR program. 

Options : The user must select the Contour option prior to plotting a view 

using Namelist SXNHJT parameters,, 

Tho options available are: 

o To punch transformed grid coordinate data on file 10. 

0 To punch Grid data on file 10 and suppress plotting. 
o v 0 punch element connection data on file 11. 

0 To punch element connection data on file 11 and suppress 

plotting. * 

° To specify alternate file numbers for output of grid and 
element data for CONTUR. 

Procedure : The user sets parameter CONTUR, = 1 or 9 at the time Namelist $ INPUT 
is specified. 

(Option 1 does not suppress the screen plot. Option 9 suppresses the 
screen plot). 

° A set of transformed grid coordinates and element connection 
numbers are written on £ilet> 10 and 11 respectively. The 
. transformed grid coordinates are the xy screen coordinates 

©£ the projected view. 

° If the output is desired in the form of punched cards the 
parameters At and IQ are set equal to the punch symbiont 
* file number (normally file I on standard UN1VAC systems). 

Limitations : When CONTUR - 1 is used 0 the multiple frame plot parameter (NPLOT 

1} is set tc terminate the run after the nest plot. IS contour 
data is required for several views s a new execution must be initiated 
for each view. 


ORIGINAL) ^ GE J® 
0 % EOOR QUALITY, 
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3.5 ALTERNATE PRINT FILE 


( Purpose : To provide a permanent record of the plotting session* 

Options : The following output is written on file 21* 

* Maximum number of elements 

* Title card 

* Namelist parameters 

* Grid and connector listing 

* Transformed grid coordinates (from local to basic coordinate 
system). 

* Number of grid points and elements in plot, 

* X,Y,Z coordinate max and min values. 

* Scale factor 

* Error messages 

Procedure : File 21 is always written. To obtain a printout) the user must 

catalog file 21 prior to execution or copy file 21 into a cataloged 
file. The file is then SYM'ed to a system printer. 

Limitations : dot Applicable* 

( 


( 
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APPENDIX A 


( 


Subroutine Descriptions 


( 


-•V 


ABSA 

Purpose: Calculates absolute value of vector VK 

External References: SQRT 

Called by: TPORM 

Calling Argument: (VK.AABS) 

VK a Vector input 

AAS3 a Absolute value of vector 

ACROSB 

Purpose: Calculates vector cross product AxB a C 

External References: None 

Called by: TFORM PAGINAL PAoli I s 

Calling Argument: (A,B, Cross) Qg j>QQR ftUAUTX 

AXIS 

Purpose: Adds three dummy bar elements for use In drawing x,y,z references 

axes* 

External Reference: None 

Called by: COMBS 

Calling Argument: (IP0INT,G1,X1,X2,X3,XMAX,YMAX)ZMAX I N0EL,N0GD I 1DG,ELTYP,IDEL) 

AXESLC 

Purpose: Calculate the co<‘-' ,< .nate position of three grid points used to 

v draw the reference axes. 

External References: None 

Called by: MAIN l 

Calling Argument: (X1,X2,X3,XMAX,YMAX,ZMAX,XMIN,YM1N,ZMIN,N0GD) 


AXIS1 


Purpose: To change max/min values for axis prepositioning (interchanging 

axes or creating a left hand coordinate system) .• 

External References: None 

Called by: MAIN1 

Colling Argument: (XAXIS.YAXIS.ZAXIS.YMAX.XMIN.YMAX.YMIN.ZMAX.ZMIN) 


AX1S2 


( 


Purpose: To change grid coordinates for axis prepositioning. 

External References: None 
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Called by: MAIN1 

Calling Argument: Xl,X2,X3,Il,l2 t Jl,XAXIS t YAXIS,ZAXIS,NGRID) 

( CINGEN 

Purpose: Read file 29 which contains thermal network data created 

by the CINGEN Program. 

External References: None 

Called by: MAIN1 

Calling Argument: (NL,NS,ICONDS,NA,NB,XHA,YNA,ZNA,XNB,YNB,ZNB) 

COMBS 

° Purpose: Counts grid points element and coordinate systems, transforms 

grid points to basic coordinate system, determines max and 
min x,y,z coordinates and renumbers grid 10*8 with internal ID*s. 
External References: AXES , TFOKM, MAXIMIM 

Calling Argument: (NOGD.LDS.Xl.XZ.XS.XMIlI.XMAX.YMIN.YMAX.ZMINpZMAX, 

NOFL.FLTYPE, IDEL , G1 , 1FOINT, CP ,NTEST) 

COUNT 

Purpose: To count the number of elements, grid points and element connection 

points to determine the minimum core size required for the given 
/ data. 

' Called by: MAIN 

Calling Argument: (NUMEL, NUMGRD.jf SIZE, SIZE, $) 

DECODE 

Purpose: Decode NASTRAN free field variables 

External References: None 

Called by: READB 

Calling Argument: (M,NUM3) 

DSCRI 

Purpose: Writes a table describing the Namelist parameter options on 

the CRT screen. 

External References: None 

Called by: MAIN1 

Calling Argument: None * 


o 


V. 


Dummy 


Purpose: To load blank common after Editor is completed 

Called by: MAIN,MAIN1 


EDITOR 

The Editor has been removed from Level 4.0. 
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GENERT (UNIVAC only) 

Putpose: To generate intermediate points for CIS3D20 so that a smooth 

curve can be fitted through the input grid points. 

External References: MAT1HV 

Called by: PL0T1 

Calling Argument: ( B , INTER V , I SKIP , 1 SIDES ) 

LOAD 

Purpose: To load values of thermal conductors to arrays for plotting. 

External References: None 

Called by: MAIN1 

Calling Argument: {NGR1D,I1,S1,NL,1N,IC0NDS,NA,NB,XNA,YNA,ZNA,XN3,YNB, 

ZNB , GTEMP ,X1 , X2 , X3) 


MAIN 


Purpose: To apportion cosnaon. into ten segments for plotting, calculates 

maximum number of elements which may be plotted, identifies file 
usage, plot title page, and clears plot buffers prior to terminat- 
ing execution. 

External References: DUMMY ,MAIN1 , COUNT , SIZ , MCORE, LCORE 

Comments: The main program vas designed to open the core size to fit 

the problem being plotted. 


MAIN1 


Purpose: To set the default values for Namelist ($INPUT) and control the 

flow of the program. 

External References: DSCRI,SPAR,RDB,TP0R,PVW,AXESLG,AXIS1,VW1,AXIS2,VW2, 

CINGEN, LOAD ,XY SAVE, EDITOR, DUMMY, STRUDL 

Called by: MAIN 

Calling Argument: <IDG,X1,X2,X3,IDEL,ELTYP,G1,1P0INT,JPV,YY) 

Comments: The communication line transmission rate can be changed through 

Namelist $INPUT prior to the first plot* See BAUD « Rate/ 10. 
The default value is 120 • 1200 Baud. Once it has been changed 
it will remain changed for the current execution. 


MATINV 


Purpose: Invert matrix for subroutine GENERT. 

External References: None 

Called by: GENERT 

Calling Argument: (A. 




8-23 


MAXS-HN 


Purpose: To find largest and smallest numbers in an array of grid coordinates X 

External References : None 

Called by: COMBS 

Calling Argument: (JJ , XX, MAX, MIN) 

MPS 

Purpose: Map symbolic for segmenting the program. 

NUMBER 

Purpose: To write a number at the grid point and/or element centroid. 

External References: None 

Called by: PLOTVG,PLOT3G,PLOT1 

Calling Argument: (XN,X,Y) 

PITTAS C 

Purpose: To switch from vector generation to hardware character generation 

after each. plot. 

External References: None 

Called by: MAXN1 

Calling Argument: None 

PLOT 

Purpose: To convert drum plotter calls to Tektronix plot calls. 

External References: None 

Called by: PLOTVG , PL0T3G f PL0T1 # YSHBOL , VIEW1 

Calling Argument: (X,Y,IP) 

PLOTVG 

Purpose: To draw a line between grid points and place a symbol and grid 

number at each grid point if requested. Used for 20 grid elements 
only. 

External References: NUMBER, SYMBOL* PLOT 

Called by: PL0T1 

Calling Argument: (C,GN0,SYMB,NL,J1,AH,AH2,PIBI2,JPT,N0,INTERV,CHAR, 

ISKIP) 

PLOT1 

Purpose: To plot finite elements and conductors and write element, grid, 

conductor, capacitor numbers at grid point and element centroids 
and to draw a symbol at the grid point if requested. 

External References: NUMBER, SYMBOL, PLOT ,XY2LAB,FL0T3G,GENERT, PLOTVG 

Called by: VIEW2 

Calling Argument: (C,GN0,EN0,KL,K1,NGRID,NND,SYMB, JK,NOEL,IBR) 


8-24 [ 3-27 


PL0T3G 


Purpose: To draw a line between grid points, and place a symbol and grid 

number at each grid point if requested. Used for all elements 
except the 20 grid isoparametric element. 

External References: NUMBER, SYMBOL, PLOT 

Called by: PL0T1 

Calling Argument*: (C S GNO,SYISB,NL,J1,AH,AH2,PIBI2,JPT,NO,CHAR) 

PVIEW 

Purpose: To read Namelist $PVIEU and calculate data for the partial view. 

External References: SEARCH 

Called by: PW 

Calling Argument: (NOEL.IDEL.Gl.NOGD.Xl.Xa.XS.IPOINT.XtaN.XMAX.YMIN, 

YHAX, ZMIN , ZMAX , YY , IDG , JP V ,N JP V , LENGB , LP 9 ) 

FVW 

Purpose: Dummy routine for calling PVIEW. 

External Reference: PVIEW 

Called by: MINI 

RDB 

Purpose; Dummy routine for calling READB 
External Refer encea: READB 

Called by : MINI 

READB 

Purpose: To read finite element data from file 4, decodes free field data, 

stores in arrays for plotting and calculates maximum and minimum 
coordinate values. 

External References; BCODE 
Called by: RDB 

Calling Argument: CNOGD 0 IDG,X1,X2,X3,XMIN,XMAX,YMIN,YMAX,2MIN,ZMAX,NOEL, 

ELTYP , IDEL , G1 , IPOINT, NCHBDY , NTEST , CP .PROGRAM, LP9 ,LENGB) 

SEARCH 

Purpose: Not fully implemented at this time. 

SPAR (UNIVAC only) 

Purpose: To convert SPAR print file (9) to NASTRAN Bulk Data format (file 4) 

External References: None 

Called by: MAIN1 

Calling Argument: None 


SYMBOL 


Purpose: To plot symbols at grid point for structural plots or at element 

centroid for thermal network plots. 

External References: PLOT 

Called by: PLOTVG , PLOTl ,PL0T3G 

Calling Argument: (X,Y,SYM) 

TFOR 

Purpose: Dummy routine for calling TFORM 

External References: COMBS 

Called by: MAIN1 

TFORM 

Purpose: To transform grid point coordinates from 3 oca! coordinate systems 

to a basic rectangular coordinate system. 

External References: ABSA,ACROSB 

Called by: COMBS 

Calling Argument: (NOGD,NCORD,X1»X2,X3,CTYPE,CP,CID,RID,A) 

Comment: Only NASTRAN C0RD2R, C0RD2C and C0RD2S coordinate definition cards 

are currently permitted. 

V1EW1 

Purpose: To calculate plot parameter, such as, center of plot and scaling 

factors. 

External References: PLOT 

Called by: VW1 

Calling Argument: (ANGLD , EL E V , ROTX , XMAX , XMIN , YMAX , YMIN , ZMAX , ZMIN , XMEAN , 

YMEAN,ZMEAN, ANGLE ,T,ROTR, Cl) 


VIEH2 


Purpose: To rotate and scale each set of grid coordinates for each element 

and call plot subroutine. 

External References: PLOTl 

Called by: VW2 

Calling Argument: (C,NGRID, XMEAN, YMEAN,ZMEAN, ANGLE, T,R0TR,NL,J1,NTEST,EN0, 

GNO , SYMB , NND , Cl , JK ,NOEL , IBR) 


VW1 


Purpose: Dummy routine for calling VIEW1 

External References: VIEW1 

Called by: MAIN1 


VW2 

Purpose; Dummy routine for calling VIEW2, 

( External References: VIEW2 

Called by: MAIN1 

XYSAVE 

Purpose: To write grid point number and X # Y location on CONTUR data file 

KC (KC - 10 9 default). 

External References: None 

Colled by: MAIN1 

Calling Argument: (KC,LAST»XDG,NGRlD,GlyC) 

XYZLAB 

Purpose: To label the trial coordinate symbol with X,Y and Z labels. 

External References: None 

Called by: PL0T1 

Calling Argument: (C1,C2,N0) 

STRUDL (UNIVAC only) 

Purpose: To read STRUDb finite element data from file 4» decodes data, 

stores in arrays for plotting and calculates max and min 
coordinate values. 

^ External References: DCODE 

Called by: MAIN1 

Calling Argument: (B(LP1) ,B(LP2) ,BaP3),B(LP4) ,B<LP5) ,B(LP7)yB(LP8)) 


